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STATEMENT OF GOVERNMENT INTEREST 
[0001] This invention was made with United States Government support from 

the Office of Naval Research (ONR) under contract number NOOO 14-02- J- 1822 and 
from the National Aeronautics and Space Administration (NASA) under contract 
numbers NCC5-563, NAG5-11819, and NAG5-12089. The United States 
Government has certain rights in the invention. 



CROSS REFERENCE TO RELATED APPLICATIONS 
[0002] The present application claims priority to U.S. Provisional Application 

No. 60/439,391 filed January 10, 2003 entitled REAL-TIME SOFTWARE 
RECEIVER which is incorporated herein in its entirety by reference. 

BACKGROUND OF THE INVENTION 
[0003] This invention relates generally to software radio receivers, and more 

specifically to a software receiver for positioning systems. 
[0004] A typical positioning system receiver, such as is used in the Global 

Positioning System (GPS), includes an antenna, a radio frequency (RF) section, a 
correlator, a signal tracking and demodulation component, and a component to 
compute the navigation solution. The antenna, which is possibly followed by a pre- 
amplifier, receives L-band GPS signals. The RF section filters and down converts the 
GHz GPS signal to an intermediate frequency in the MHz range. The RF section also 
digitizes the signal. The correlator separates the down-converted signal into different 
channels (ten or more in modem receivers) allocated to each satellite. For each 
satellite, the correlator mixes the Doppler-shifted intermediate frequency signal to 
baseband by correlating it with a local copy of the carrier replica signal and it 
distinguishes the particular satellite by correlating the signal with a pseudo-random 
number (PRN) code. Software routines cause the carrier replica and PRN replica 
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signals to track the actual received signal, extract the navigation message, and 
compute the navigation solution. 

[0005] Baseband mixing is a multiplication of an input signal by a complex 

exponential where the frequency of the complex exponential approximately matches 
that of the input signal The resultant signal is centered at baseband. A complex 
signal can be broken down into cosine and sine signal components, resulting in 
separate in-phase and quadrature components. The frequency of the baseband mixed 
signal must be controllable to within a few millihertz in the case of a phase-locked 
loop for use in a precision navigation system, and the baseband mixed signal must 
have a continuously varying phase. In a hardware correlator, local oscillators generate 
cosine and sine signal components that have precise frequency control and a 
continuous phase. Generating cosine and sine signal components on the fly with the 
correct frequency and phase is too time consuming to be feasible for a software 
correlator. Instead, the software correlator generates cosine and sine signal 
components on a grid of frequencies off-line. These signal components must be stored 
on a time grid of points sampled at the RF front-end sampling frequency, for example, 
at 5.714 MHz for one particular RF front-end hardware configuration, and the signals 
must last for a typical accumulation period , e.g., for a 0.001 second coarse/acquisition 
(C/A) PRN code period when working with GPS LI civilian signals. It takes tens of 
gigabytes of memory or more in order to brute-force store all frequencies on a one 
mHz grid ranging from -10 KHz to +10 KHz, which is the needed frequency range 
when tracking GPS satellites from a terrestrial receiver, and additional storage is 
required to store a grid of possible starting phases at each frequency point. 
[0006] PRN code mixing is a muhiplication of a baseband mixed signal by a 

prompt +1/-1 PRN code or by a +2/0/-2 early-minus-late PRN code, where the code 
timing and frequency approximately match that of the input signal. The resultant 
signal is a constant in the case of prompt PRN code mixing, and an approximately 
linear function of the code timing error in the case of early-minus-late mixing. A 
receiver accumulates both of these correlation outputs. The magnitude of the prompt 
accumulation indicates signal strength and whether a signal has been detected, and its 
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in-phase (real) and quadrature (imaginary) components are used to measure carrier 
phase and Doppler shift. The magnitude of the early-minus-late accumulation 
measures the code timing error; it will be zero when the timing error is zero. 
[0007] The code phase of the baseband mixing signal must be controllable to 

within a percent or less of a PRN code chip for use in a precision navigation system. 
In a hardware correlator, local oscillators generate the prompt and early-minus-late 
PRN code replicas. A software correlator can either compute and store PRN code 
replicas, or compute them in real-time. 

[0008] The current Global Positioning System is slated to realize expanded 

capabilities that include new civilian codes on the L2 frequency, a new L5 frequency, 
and new codes (M-code, CL and CM codes) on the L2 frequency. Some of these 
upgrades are slated to start within one to three years. A hardware correlator requires 
hardware modifications in order to use these new signals. In the near term, a receiver 
designer will be faced with a complex trade-oflf in order to decide whether the extra 
complexity is worth the improved performance that will accrue only very slowly as 
new GPS satellites replace older models. One way to avoid the complex trade-off is to 
use a software receiver that can receive and process new signals without the need for a 
new correlator chip set. 

[0009] A software receiver is flexible because its software components can be 

easily modified. One application of a software receiver is to merge together numerous 
devices that use Mdreless digital communication protocols to form a single device. For 
example, a cell phone, GPS receiver, and Personal Data Assistant (PDA) could 
become a single device that plays the role of all three. Another use of a software 
receiver is to shorten development and to-market times for new wireless devices. For 
example, as new frequencies and codes are added to GPS, a software receiver having a 
software correlator simply needs to be reprogrammed, while a hardware approach 
would require a brand new correlator chip design. New PRN codes can be used simply 
by making software changes. Thus, software receiver technology lessens the risks 
involved for designers during the period of transition to the new signals. Furthermore, 
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a software receiver could be reprogrammed to use the Galileo system (European GPS) 
or GLONASS (Russian GPS). 

[00010] In the recent past, GPS software receivers have been developed that 
either post-process stored signals or operate in real-time. Previous real-time software 
receivers fiinction with a limited number of channels (4-6) or require high-end 
computer speeds or digital signal processor (DSP) chips such as are disclosed in Real- 
Time GPS Software Radio Receiver, Akos et al., ION NTM 2001, 22-24 January 2001 , 
Long Beach, CA, pp. 809-816 {Akos 2001a), and Global Positioning System Software 
Receiver (gpSrx) Implementation in Low Cost/Power Programmable Processors, 
Akos et al., ION GPS 2001, 11-14 September 2001, Salt Lake City, UT, pp. 2851- 
2858, both incorporated herein in their entireties by reference. 
[00011] Therefore, it is an object of the present invention to create a software 
receiver that operates in real-time and is not restricted to a severely Hmited number of 
channels or to a very fast processor. 

[00012] Another object of the present invention is to minimize the number of 
sine and cosine signal components that must be stored. 

[00013] A fiirther object of the present invention is to process incoming signals 
through bit-wise parallelism. 

[00014] A still fiirther object of the present invention is to process over-sampled 
signals by use of bit-wise parallelism. 

[00015] A still further object of the present invention is to use very long over- 
sampled PRN codes efficiently in a bit-wise parallel software receiver. 



SUMMARY OF THE INVENTION 
[00016] The objects set forth above as well as fiirther and other objects are 
addressed by the present invention. The solutions and advantages of the present 
invention are achieved by the illustrative embodiment described herein below. 
[00017] The software receiver system and method of the present invention 
enable the efficient execution of a set of algorithms, that perform software correlation 
on data sampled fi-om incoming channels, on a general purpose processor. The system 



5 



and method of the present invention provide for either PRN code storage or 
computation of PRN codes in real-time. PRN code storage is appropriate for PRN 
codes that have short periods, such as the GPS coarse/acquisition codes, which are 
1023 chips long. In this case, the system and method of the present invention pre- 
compute over-sampled replicas of entire PRN code periods and store them for orderly 
and efficient retrieval, such as in a table. This table can include a selection of code 
start times as measured relative to the sample times at which RF data are available 
from the receiver front end. There is a separate table for each unique PRN code. 
[00018] The system and method of the present invention can also generate over- 
sampled versions of the prompt and early-minus-late PRN codes in real-time through 
use of an over-sampling fiinction described herein. The values of the over-sampling 
function can be located in a specially designed table that can be generic across PRN 
codes. The length of the specially designed table can be independent of the length of 
the PRN code whose replica is being used to process a given received signal. The 
system and method of the present invention include techniques for efficiently 
calculating indices into the specially designed table that enable rapid, real-time table 
look-up. 

[00019] The system of the present invention includes a software correlator that 
can mix the received signal to baseband, compute baseband/PRN correlations through 
bit-wise parallelism and look-up tables using either the tabulated or real-time- 
generated PRN codes, and compute accumulations through bit-wise parallelism and 
processor instructions or look-up tables. Bit-wise parallelism allows the processing of 
multiple data samples simultaneously as the multiple bits of a given word of computer 
data. For example, for 32-bit words, the software correlator can process up to 32 
samples at a time. Bit-wise parallelism can optimally operate when each signal in 
question can be represented by only a few bits, which is normally the case in RF 
digital signal processing of navigation signals. 

[00020] The bit-wise parallel operations of the present invention can save 
computation time in comparison to integer mathematical correlation operations. If, for 
example, four accumulations are required per sample, integer mathematics requires six 
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multiplications and four additions per sample (except for the last sample). At a 
sampling rate of, for example, 5.714 MHz this translates into 57,140 integer 
operations per PRN code period. In the illustrative embodiment, 33,500 bit-wise 
parallel operations are necessary per PRN code period when the RF signal has a 2-bit 
representation. This operation count is further reduced to approximately 16,750 bit- 
wise parallel operations per PRN code period when the RF signal has a 1-bit 
representation. Thus, there can be a savings of almost a factor of two to almost a 
factor of four in the operation count. 

[00021] The system and method of the present invention also include a table of 
pre-computed baseband mixing sine waves, algorithms that can produce correlation 
accumulation outputs that are equivalent to what would be produced by a continuously 
variable sine wave, and a method of use of the table and algorithms. Thus, in the 
present invention, a relatively small set of sine wave values need to be pre-computed 
and saved, which can conserve computer memory and processing time. 
[00022] The present invention also includes a system and method for tracking 
the phase of PRN code replicas in software in order to track the timing of any given 
"chip" of the PRN code replica as measured with respect to a pre-specified set of 
sample times at which the basic raw data comes out of the RF front end (a chip is an 
element of a PRN code). The PRN code phase is kept track of via a variable for each 
channel, that indicates the PRN code start time with respect to the RF sample times. 
The system and method of the present invention allow for the synchronization of the 
measurements of PRN code phase, carrier phase, and carrier frequency for each 
satellite relative to these sample times. 

[00023] The method for tracking the phase of each PRN code replica and the 
phase of each carrier replica includes the steps of latching all the C/A code phases, 
carrier phases, epoch counters, and carrier frequencies for each satellite at a pre- 
specified time, and computing the pseudo range to each satellite using the C/A code 
phase and epoch counters. The method also includes the step of tracking and updating 
code and carrier phases by estimating code chipping rate and carrier Doppler shift 
inputs. The method fiirther includes the step of computing the code phase at the pre- 
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specified time for each satellite as a function of the updated code chipping rate and the 
pre-specified time. The method fiirther includes the step of computing the carrier 
phase at the pre-specified time as a function of the updated carrier phase, the Doppler 
shift, and the pre-specified time. The timing of the PRN code phase (or chip location) 
is the most fundamental of GPS measurements for use in navigation data processing. 
The monitoring of these times in software allows complete control of the precision 
with which they can be measured, and it allows precise synchronization of these times 
with the measurement times of data from other sensors, such as inertial measurement 
units. This feature gives an enhanced ability to develop what are known as deeply 
coupled systems that must fiise GPS data with data from other types of sensor 
systems. 

[00024] The software correlator of the present invention can advantageously be 
easily adapted to accept signals at any frequency, new PRN codes, or even signals for 
different types of devices. Thus, the same processing hardware could use the software 
correlator to implement such devices as a GPS receiver, a cell phone, or both. To 
allow for new codes, new frequencies, and new types of functionality, small changes 
can be made in the software correlator, or different versions of the software correlator 
can be run on the same processor. Hardware-correlator-based receivers of the prior art 
can deal only with frequencies and PRN codes that are hard-wired into their designs. 
Also, the system and method of the present invention could be implemented within 
systems such as GLONASS receivers, cell phones and cell base stations, pagers, 
wireless Ethernet (e.g. 802.11x standards), Bluetooth™, Blackberry® wireless internet 
devices, and satellite radio/phones (e.g. INMARSAT®). In fact, the system and 
method of the present invention are applicable to any sort of telecommunication 
system/device that uses spread spectrum, code division multiple access (CDMA) PRN 
codes for the transmission of information, either wired or wireless. 
[00025] For a better understanding of the present invention, together with other 
and further objects thereof, reference is made to the accompanying drawings and 
detailed description. The scope of the present invention is pointed out in the appended 
claims. 
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DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 
[00026] FIG. 1 is a schematic block diagram of the hardware environment of a 
typical software receiver; 

[00027] FIGs. 2A and 2B are schematic diagrams of bit-wise mappings of 

signal and carrier replica sign and magnitude bits to computer data words; 

[00028] FIG. 2C is a graphic representation of a plot of bit- wise parallel radio 

frequency signal and PRN code replica storage and mixing; 

[00029] FIG. 2D is a graphic representation of a plot of sections of prompt, 

early, late, and early-minus-late PRN code signals and 16-bit word representations of 

their over-sampled equivalents; 

[00030] FIGs. 3 A and 3B are data flow diagrams illustrating the bit-wise 
parallelism process (replicated twice, once for the in-phase carrier replica and once for 
the quadrature carrier replica) of the present invention; 

[00031] FIG. 3C is a graphic representation of a plot of a prior art optimal 2-bit 
representation of a sine wave presented to enhance the reader's understanding of the 
present invention; 

[00032] FIGs. 4A and 4B are flowcharts of the method for computing 
correlation accumulations through bit-wise parallel computations of the present 
invention; 

[00033] FIG. 5 is a schematic diagram of a look-up table layout as a function of 
code time offset and chip bit pattern; 

[00034] FIG. 6 is a graphic representation of a plot illustrating the timing 
relationship between data sample words and the sequence of prompt code chips that 
defines an accumulation interval; 

[00035] FIG. 7 is a flowchart of the method for computing bit-wise parallel 
representations of the over-sampled prompt PRN code replica and the over-sampled 
early-minus-late PRN code replica for an entire accumulation interval using the real- 
time over-sampled PRN code generation algorithm. 
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[00036] FIG. 8 is a graphic representation of a plot that illustrates the location 
in time at which the code phase of each signal is computed; and 
[00037] FIGs. 9A and 9B are graphic representations of plots of correlations of 
the true sampled code with prompt (FIG. 9A) and early-minus-late (FIG. 9B) versions 
of the true and table look-up codes, the latter being generated by the new real-time 
over-sampled PRN code generator. 

DETAILED DESCRIPTION OF THE INVENTION 
[00038] The present invention is now described more fiilly hereinafter with 
reference to the accompanying drawings, in which the illustrative embodiment of the 
present invention is shown. The following configuration description is presented for 
illustrative purposes only. Any computer configuration satisfying the speed and 
interface requirements herein described may be suitable for implementing the system 
of the present invention. The equations herein are stated in general terms, but have 
parameters that are specific to the GPS LI C/A signal for illustrative purposes only. 
For example, the 0.001 sec. accumulation interval seen in many of the equations is the 
nominal C/A code period. Also, the C/A PRN code of the illustrative embodiment can 
be replaced by the PRN code of any other CDMA signaling system. 
[00039] By way of introductory explanation, RF signal processing equations 
and terms are herein provided. The time-domain LI C/A signal received from, for 
example, a satellite, is represented by: 

where is the sample time, Aj is the amplitude, Dj^^ is the navigation data bit, Cj[t] 
is the C/A code, r^^. and r^j^^, are the start times of the received and ^ + r' C/A 
code periods, cDj^ is the intermediate frequency corresponding to the LI carrier 
frequency, ^j(t^) is the carrier phase perturbation due to accumulated delta range, 
is the receiver noise, and the subscript j refers to a particular GPS satellite. The 
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summation is over all visible GPS satellites. The negative sign in front of ^(/, ) comes 

from the high-side mixing that occurs in the RF front-end that has been used in the 
illustrative embodiment. The signal in equation (1) is the output of a typical RF front 
end. 

[00040] A GPS receiver works with correlations between the received signal 
and a replica of it. The correlations are used to acquire and track the signal. The 
replica is composed of two parts, the carrier replica and the C/A PRN code replica. 
Two carrier replica signals are used, an in-phase signal and a quadrature signal. When 
mixed with the received signal from the RF front end they form the in-phase and 
quadrature baseband mixed signals represented by: 

y&iO = -C,[0.001 (/'"^^* ) ]sm^^t, -[^.+Soop„tit, -f } (3) 

Tjk^l'-'Tjk 

where equations (2) and (3) apply during the A:* C/A code period. In these equations 
Tjk and Tjk+\ are the receiver's estimates of the start times of the k^^ and code 
periods, ^.^ is the estimated carrier phase at time zjk , and d>Domk *s the estimated 
carrier Doppler shift during the k'^ code period. 

[00041] A typical receiver computes the estimates Zjk > Tjk+i , (^^^ , and Soomk 

by various conventional means that are described in GPS Receivers, A. J. Van 
Dierendonck, Global Positioning System: Theory and Applications. B. W. Parkinson 
and J. J. Spilker, Jr., Eds,, vol. I, American Institute of Aeronautics and Astronautics, 
1996, Chapter 8, pp. 329-406 {Dierendonck\ incorporated herein in its entirety by 
reference. These include open-loop acquisition methods and closed-loop signal 
tracking methods such as a delay-locked loop to compute fy* and ijk+x and a phase- 
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locked loop or a frequency-locked loop to compute ^.^ and a>Domk The software 

receiver developed herein uses conventional techniques for forming these estimates. 
[00042] Both prompt and early-minus-late correlations are needed to track the 
carrier frequency, carrier phase, and code phase in a GPS receiver. A typical receiver 
uses the PRN code and carrier replicas to compute the following in-phase and 
quadrature correlation accumulations: 



i=ii. T jk+\ T Jk 



h'^^k ^ _|_ A — '* 

QA^) = - Z X^,)C,[0.001( \ ^/ )] sin^^t, -[^^^ + a>DoppAt, - tj,)] } (5) 

1=/,. T jk+\ T jk 

where ik is the index of the first RF front-end sample time that obeys ij^ ^ and 

Nk-^\ is the total number of samples that obey i j]^ < r,- < ryj^+i . The time offset A 

causes the replica PRN code to play back early if it is positive and late if ^ is 
negative. The prompt correlations are defined by equations (4) and (5) v^th A = o. 

The early-minus-late correlations are ^^^^^^"^^^^ ^em/^^) 

Qjk (^*m//2) - Qj^ ^ ^j^^ spacing between the early and late PRN 

carrier replicas. The present invention described herein is an efficient technique for 

the receiver to accumulate and in software. 

[00043] Referring now to FIG. 1, the operational platform of the software 
receiver 10 of the present invention includes an antenna 1 1, conventional RF front-end 
13, a data acquisition (DAQ) system 17, a microprocessor 16, a software correlator 19, 
and application-specific code 15. Conventional RF front-end 13 interfaces v^th 
antenna 1 1 and with (DAQ) system 17. DAQ system 17 includes a system of shift 
registers and a data buffer. Microprocessor 16 executes software correlator 19, which 
includes a set of specially developed bit-wise parallel algorithms, and application- 
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specific code 15, such as the GPS navigation and tracking fiinctions. In the 
illustrative embodiment, conventional GPS software functions (signal tracking, data 
extraction, navigation solution, etc.) are provided by the MITEL® GPSArchitect 
software ported to RTLINUX® (see A Coming of Age for GPS: A RTLINUX BASED 
GPS RECEIVER, Ledvina et al.. Proceedings of the Workshop on Real Time 
Operating Systems and Applications and Second Real Time Linux Workshop (in 
conjunction with IEEE RTSS 2000) November 27-28, 2000), but can be provided by 
any equivalent configuration. 

[00044] Continuing to refer to FIG, 1, conventional RF front-end 13 can, for 
example, be a MITEL® GP2015 RF front-end, which down converts the nominal 
L57542GHZ GPS signal 12 to an intermediate frequency of (88.54/63) x 10^ Hz 
1.4053968254 MHz and then performs analog-to-digital conversion. The resultant, 
digitized signal data 21 has a pre-determined number of bits/sample, such as two 
binary bits/sample, a sign bit and a magnitude bit, or one bit/sample. The shift 
registers in the DAQ system 17 parallelize the magnitude and sign data bit streams 
into separate words, which the DAQ system 17 reads into the memory of 
microprocessor 16 using DMA. To make the process of reading data into the 
microprocessor 16 more efficient and to prepare for efficient correlation calculations, 
DAQ system 17 can read a pre-specified number of bits of buffered samples, such as 
thirty-two bits, at a time. The exemplary thirty-two bits include sixteen sign bits and 
sixteen magnitude bits. 

[00045] Referring now to FIGs. 1 and 2A, the shift registers in DAQ system 17 
(FIG. 1) buffer signal data 21 (FIGs. 1 and 2A) and pack signal sign 21 A (FIG. 2A) 
and signal magnitude 2 IB (FIG. 2 A) into separate words, that represent the integer 
values ±1 and ±3 as is shown in Table 1 . In the case of a 1-bit signal, the bit stream 
representing the samples is packed into successive words to prepare the signal for bit- 
wise parallel processing. DAQ system 17 also provides for accurate timing by, for 
example, synchronizing signal sign 21 A and signal magnitude 2 IB to a (40/7) x 10^ 
Hz ^ 5.714 MHz clock signal, which can be, for example, a third output from 
conventional RF front-end 13 (FIG. 1). DAQ system 17 can convert the 5.71424MHz 
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clock signal down to 357. 14KHz by use of, for example, a divide-by-16 counter for a 
16-bit word, which can provide a signal indicating when the buffer is full. DAQ 
system 17 can use any method for providing a buffer full indication. 



Signal Sign 21A 


Signal Magnitude 21B 


RF Signal Value 


0 


0 


-1 


0 


1 


-3 


1 


0 


+1 


1 


1 


+3 



Table 1. The sign and magnitude combinations of the sample RF output of the 
conventional RF front-end and their corresponding values. 



[00046] With further reference to FIG. 1, in the illustrative embodiment, the 
DAQ system 17 can consist of an interface card and driver software that can be 
compatible with, for example, a 1.73GHz AMD ATHLON® processor running 
RTLINUX®, but could be compatible with any operating system and any processor 
that can accommodate real-time operations. The interface card can, for example, be a 
NATIONAL INSTRUMENTS® PCI-DIO-32HS digital I/O card. Pertinent features 
of this card are the thirty-two digital input lines, DMA, and availability of a driver for 
RTLINUX®, perhaps gotten from the suite of open source drivers and application 
interface software for interface cards known as COMEDI (COntrol and MEasurement 
and Device Interface). Modifications to the conventional COMEDI driver for the PCI- 
DIO-32HS card include increasing the number of input bits from sixteen to thirty-two, 
enabling DMA, and modifying the driver to support continuous interrupt-driven 
acquisition. 

[00047] With still further reference to FIG. 1, microprocessor 16 can be, for 
example, a 1.73 GHz AMD ATHLON™ processor running the RTLINUX® operating 
system, but any operating system and processor that can accommodate real-time 
operations can be used. Low latency interrupt responsiveness, the ability to execute 
threads at regular intervals, with the kernel having a possibility of being the lowest 
priority thread, and reUable execution of time-critical code are among features of an 
operating system that could enhance the performance of the system of the present 
invention. The use of RTLINUX® is presented herein for illustrative purposes only. 
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[00048] Continuing to refer to FIG. 1, analogous to a hardware correlator that 
takes input directly from the RF front end in serial fashion, software correlator 19 
reads from a shared memory buffer that both software correlator 19 and DAQ system 
17 can access, the former to read data, and the latter to write data. The shared memory 
buffer can be implemented as a DMA memory space and a circular buffer. In the 
illustrative embodiment in which the system and method of the present invention are 
used in a GPS (or similar) environment, microprocessor 16 can store the most recent 
twenty-one milliseconds of signal data 21 (FIGs. 1 and 2 A) in the circular buffer, but 
could store more or less. The present invention does not fix the size of the circular 
buffer, nor the amount of RF data that can be stored there. The circular buffer allows 
the processing of code periods that start and stop at different times for different 
satellites during different iterations of a regularly scheduled program thread. DMA 
memory space can be written to directly by DAQ system 17 using a DAQ software 
driver, which fills the circular buffer. Communication between software correlator 1 9 
and application-specific code IS can be performed using operating system-provided 
shared memory capability. For example, the mbuff driver, included with RTLDSfUX®, 
can be used to create and manage this shared memory space. Any memory 
management system that accommodates real-time processing can be used. If the mbuff 
driver is used, kernel modules can share memory and the kernel can be restricted from 
swapping the shared memory space to long-term storage. 

[00049] Continuing with the analogy to hardware correlation, and still referring 
primarily to FIG, 1, in hardware correlation, the correlator receives frequency and 
phase information from tracking and acquisition loops that are part of application 
code, and Numerically Controlled Oscillators (NCOs) generate signals that correspond 
to the written frequencies and phases. In contrast, software correlator 19 includes 
simulated carrier and code NCOs that receive their frequency commands from 
application-specific code 15. Software correlator 19 uses these frequency commands 
to reconstruct carrier replica signal 25 (FIG. 3 A) and prompt PRN code 29 and early- 
minus-late PRN code 35 (FIG. 3 A) which it mixes with the signal data 21 (FIG. 2 A) 
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resulting in fiilly mixed prompt integrand 31 and fiilly mixed early-minus-late 
integrand 33 (FIG. 3 A). 

[00050] To further continue the analogy, a hardware correlator generates in 
real-time a particular C/A code replica at the correct Doppler shifted frequency and 
phase. In contrast, software correlator 19 can generate C/A codes off-line and store 
them in a memory table, the pre-computed over-sampled PRN code table 28 (FIG. 
3 A). The pre-computed over-sampled PRN code table 28 is used to select PRN codes 
with the correct timing relationship to the sample times of signal data 21 (FIG. 3 A). 
The codes are then used to form correlations with baseband mixed signals 23 (FIG. 
3 A), the result from which is summed to produce the standard in-phase and 
quadrature, summed prompt accumulation 45 (FIG. 3B) and summed early-minus-late 
accumulation 47 (FIG. 3B) that are equivalent to what would be produced by a 
continuously variable sine wave. These are provided to application-specific code 15, 
such as conventional GPS software that executes signal tracking and navigation 
functions. In a second approach, software correlator 19 can generate the PRN carrier 
replicas on-line at the code chipping rate and can use tabulated functions to re-sample 
the code at the sample rate of the RF front-end for purposes of calculating 
accumulations. Real-time over-sampled PRN code generator 3 OA (FIG. 3 A) is used in 
place of pre-computed over-sampled PRN code table 28 (FIG. 3 A) in this latter 
approach. This latter method can be used with longer PRN codes, such as the new 
civilian GPS L2 CL codes. 

[00051] With still fiirther reference to FIG. 1, since the received LI raw signal 
12 can have an uncertain carrier phase, software correlator 19 computes both in-phase 
(/) and quadrature (0 accumulations, as defined in equations (4) and (5). Software 
correlator 19 begins the accumulation process by using carrier replica signal 25 (FIG. 
3 A), which it gets from pre-stored carrier replica table 30 (FIG. 3 A). The carrier 
replicas in this table fall on a rough frequency grid, and they all start with a particular 
phase, for example a phase of zero. The baseband mixing process involves selecting a 
carrier replica signal 25 (FIG. 3 A) from carrier replica table 30 (FIG. 3 A) that is at the 
frequency that is as close to "ideal" as possible. In the case of a 175 Hz grid spacing, 
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the baseband mixing process selects a signal that is maximally within =bg7.S Hz of the 
ideal signal. The rough frequency grid can have a spacing of, for example, 175 Hz but 
could be larger or smaller depending on (a) the frequency range needed to cover, for 
example, ±10 Khz, (b) the amount of space available for storing pre-computed signals, 
and (c) other design decisions. The pre-computed signals in carrier replica table 30 
(FIG. 3 A) each may occupy 1 80 32-bit words in order to be guaranteed to cover the 
full 5,714 RF front-end samples that occur in one PRN code period for any possible 
code period start time within the thirty-two samples of the initial word. Thus, 
180*4=720 bytes could be required for each bit of each pre-computed carrier replica 
signal 25 that is stored in the table. The sine and cosine waves of carrier replica 
signals 25 (FIG. 3 A) each have 2-bit representations, which translates into a storage 
requirement of 2880 bytes for the carrier replica signals 25 at a given Doppler shift. 
There are 1 15 Doppler shifts that may be stored in order to cover the -10 KHz to +10 
KHz range with a 175 Hz grid spacing. This translates into 323 Kbytes of storage for 
all of the carrier replica signals 25. This approach avoids the need to pre-compute sine 
waves with a prohibitively large number of possible frequencies and phase offsets and 
it avoids the need to compute sine waves in real-time. Instead, the errors created by 
using pre-defined sine wave replicas are compensated for by post-processing 
calculations, as described below. 

[00052] In any case, and continuing to refer to FIG. 1, the resulting 
accumulations are 

ik+^t ^ + A - X 

i=1t Tjk+\~Tjk 

i^+ATj. ^ + A 

e«*(A) = -i; MOC,[0.001(-L LE)^sin[(a,^-a>^^)it,-t,^,)] (7) 

where 0?^.^^ is the grid frequency that is closest to the estimated frequency d>Doppjk and 
where t^^^^ is the time at which this carrier replica signal 25 (FIG. 3 A) has zero carrier 
phase. Software correlator 19 rotates these accumulations in order to create accurate 
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approximations of what would have been computed had the estimated carrier phase 
time history in equations (4) and (5) been used: 

/^,(A) = /^,(A)co5(A^t«.) + e«.(A)5/«(A«*^,) (8) 

e,,(A) = -I^,{A)sm(A^„^) + Q^,iA)cosiA^^^) (9) 
where A^^^^^ is the average phase difference between the grid carrier phase and the 
estimated carrier phase averaged over the accumulation interval: 

Note that equations (8), (9), and (10) are an illustrative example of how software 
correlator (19) can rotate its / and Q accumulations in order to correct for phase and 
frequency errors in its table of pre-computed carrier replica signals. There exist other 
formulas that yield equivalent results, and this patent disclosure covers all such 
techniques. 

[00053] The validity of equations (8) and (9) is dependent on the assumption 
that 

l-COS[^{0)^, - d)Doppik){Tjk.l - Tjk)] « 1 (1 1) 

For example, a 175 Hz grid spacing and a nominal C/A PRN code period of 0.001 sec 
yields a value on the left-hand side of inequality (1 1) of 0.04, which respects the 
assumed limit. 

[00054] Note that equations (8) and (9) can be derived from equations (4) and 
(5) as follows. First, the carrier phase of the grid signal in the arguments of the cosine 
and sine terms of equations (6) and (7) are added to and subtracted from the arguments 
of the cosine and sine terms in equations (4) and (5). Next, trigonometric identities are 
used to split the resulting cosine and sine terms into sums of products of cosine and 
sine fiinctions. In each product, one of the terms involves an argument like the 
arguments in the trigonometric terms in equations (6) and (7). The other trigonometric 
terms are then approximated by either cos{A^^^f^) or 5/«(A^^^^) . These 

approximations are valid because of the inequality in equation (1 1) and because the 
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average of sin{{a)gjk G}Doppjk)[h '^(^Jk +^yVk+l)]} over the accumulation interval 
is zero. 

[00055] A decrease in the carrier to noise ratio C/No, which characterizes the 
receiver's sensitivity, is caused by the use of an inexact baseband mixing frequency. 
The worst-case decrease is expressed as a fiinction of the frequency grid spacing A/* 
and is given by 

ASNR = 20feg„(^!!^) (12) 

where A/' is in units of Hz, and Tis the integration period. Thus, a A/* of 175 Hz 
causes a worst-case C/No loss of 0. 11 dB for T = ft 001 sec. 
[00056] Referring now to FIGs. 2A, 2B, and 3 A, PRN codes (composed of 
prompt PRN codes 29 (FIG. 3A) and early-minus-late PRN codes 35 (FIG. 3A)) are 
either pre-computed or generated in real-time. Pre-computing involves, for each 
satellite, computing an entire PRN code, storing the PRN code appropriately for easy 
retrieval, and referencing the PRN code, possibly by means of indices that are 
computed based on, for example, the incoming RF signal data 21 (FIGs. 2A and 3 A). 
Pre-computing can be most advantageously used when the PRN code is not very long. 
Generating PRN codes in real-time can be a more appropriate solution when the PRN 
codes are very long (and thus would require an unacceptable amount of storage), or 
perhaps when too many PRN codes are required for the amount of storage available, 
or for any other reason, but real-time PRN code generation can entail an additional 
computational cost. Both pre-computing and real-time determination of PRN codes 
are described herein with respect to a bit-wise parallel implementation. 
[00057] Continuing to refer primarily to FIGs. 2A, 2B, and 3 A, in order to 
perform bit-wise parallel operations, software correlator 19 (FIG. 1) stores pre- 
computed carrier replica sign 25 A (FIGs. 2B and 3 A) and carrier replica magnitude 
25B (FIGs. 2B and 3 A) in data words. Simple representations of signal data 21 (FIGs. 
2 A and 3 A) and carrier replica signal 25 (FIGs. 2B and 3 A) in terms of one, two, or 
more bits are suitable for using bit-wise parallelism to perform the calculations 
described herein. Bit-Avise parallel operations work with representations of the data 
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that store successive samples in successive bits of a word. For example, thirty-two 
samples (bits) of the RF front-end output are stored in two AT = 32-bit words, signal 
sign 21 A (FIGs. 2A and 3A) and signal magnitude 21B (FIGs. 2A and SA), or simply 
a single 32-bit word if signal data 21 consists of a single data bit. Carrier replica sign 
25 A (FIGs. 2B and 3 A) and carrier replica magnitude 25B (FIGs. 2B and 3 A) are 
stored, for example in tables, in separate words, with each 32-bit word storing thirty- 
two sign or magnitude bits that tabulate to thirty-two successive samples of the 
corresponding cosine or sine wave. Similarly, tables can store prompt PRN code 29 
(FIG, 3 A) and early-minus-late PRN code 35 (FIG. 3 A), which are composed of 
prompt PRN code sign 29 A (FIG. 3 A), early-minus-late PRN code sign 35 A (FIG. 
3 A), and early-minus-late PRN code zero mask 35B (FIG. 3 A). The data words that 
comprise the bit-wise parallel representations of these three signal types, the original 
RF signal data 21 (FIGs. 2A and 3 A), the carrier replica signal 25 (FIGs. 2B and FIG. 
3 A), and the de-spreading prompt PRN code 29 (FIG. 3A) and early-minus-late PRN 
code 35 (FIG. 3A), are the inputs to the calculations of software correlator 19 (FIG. 1). 
[00058] Further continuing primarily to refer to FIGs. 2A, 2B, and 3 A, many 
intermediate calculated quantities and at least three types of intermediate signals are 
also stored in bit-wise parallel format. First there are the in-phase and quadrature 
baseband mixed signals 23 (FIG. 3A), whose 3-bit representations for the illustrative 
embodiment are stored as baseband mixed sign 23 A (FIG. 3 A), baseband mixed high 
magnitude 23B (FIG. 3 A), and baseband mixed low magnitude 23C (FIG. 3A). The 
second bit-wise parallel signal type is the ftiUy mixed integrand, of which there are 
four signals: in-phase and quadrature folly mixed prompt integrand 31 (FIG. 3 A) and 
in-phase and quadrature folly mixed early-minus-late integrand 33 (FIG. 3 A). The 
former are stored as 3-bit representations in the illustrative embodiment as folly mixed 
prompt integrand sign 31A (FIG. 3 A), folly mixed prompt integrand high magnitude 
3 IB (FIG. 3 A), and folly mixed prompt integrand low magnitude 31C (FIG. 3 A). The 
latter are stored as 3.5-bit representations in the illustrative embodiment as folly mixed 
early-minus-late integrand sign 33A (FIG. 3A), folly mixed early-minus-late integrand 
high magnitude 33B (FIG. 3 A), folly mixed early-minus-late integrand low magnitude 
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33C (FIG. 3A), and fully mixed early-minus-late integrand zero mask 33D (FIG. 3A). 
This representation is called a 3.5-bit representation because the sign, high-magnitude, 
and low-magnitude bits are ignored if the corresponding zero mask bit has the value 
zero. The third bit-wise parallel signal type is a value word, of which there are two 
types: prompt integrand value words 27 (FIG. 3B) and early-minus-late integrand 
value words 37 (FIG. 3B). Each folly mixed integrand is used to construct value 
words, one word for each possible value that the integer integrand can take on. There 
are eight possible values for the integrands of the illustrative embodiment: -1, -2, -3, - 
6, 1,2, 3, and 6 for the in-phase and quadrature folly mixed prompt integrands 3 1 
(FIG. 3 A) and -2, -4, -6, -12, 2, 4, 6, and 12 for the in-phase and quadrature folly 
mixed early-minus-late integrands 33 (FIG. 3A). Each bit-wise parallel value word 
contains a one bit for each sample time when the integrand value equals the value of 
the value word, but it cont^ns a zero bit for all other sample times. The storage of 
raw data and intermediate results in bit-wise parallel format allows the EXCLUSIVE 
OR operations that are involved in mixing to operate on thirty-two samples at a time if 
microprocessor 16 (FIG. 1) has a bit-wise EXCLUSIVE OR command. Other bit- 
wise commands are used to perform additional software correlation operations in 
parallel on sets of two thirty-two samples. 

[00059] At this point, the problem of over-sampling is introduced. Referring 
now to FIGs. 2C and 2D, the problem of over-sampling is illustrated with respect to 
bit-wise parallelism as follows. There is normally more than one RF data sample per 
PRN code chip. The three successive -1 values 73 (FIG. 2C) at sample times ti to all 
occur during the same PRN code chip as do the four successive +1 values 75 (FIG. 
2C) at times t4 through Z?- The difference in the number of samples for the two code 
chips arises because the PRN code chip period is not an integer multiple of the sample 
period. Analogously, referring to FIG. 2D, where sample interval Ats 63 is less than 
actual PRN code chip length /kc 65, over-sampling is indicated because the RF 
sampling frequency fs = is greater than the PRN code chipping frequency fc = 
}/Atc, PRN codes for CDMA signaling are sequences of +1 and -1 values, the 
elements of which are chips. Over the time intervals of interest, a carrier replica 
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progresses through its chips at a constant chipping rate of/? = l/Atc chips/second. The 
time interval ^d/e is the actual PRN code chip length 65 (FIG. 2D). Software correlator 
19 (FIG. 1) normally receives PRN code, and attempts to align it with the prompt 
replica version of the code, prompt PRN code 29 (FIG. 2D). It makes use of the 
signal's correlation with prompt PRN code 29 (FIG. 2D) and with early-minus-late 
PRN code 35 (FIG. 2D) in order to determine a chipping rate fc that tends to align 
prompt PRN code 29 (FIG. 2D) as desired. Conventional methods for determining/: 
are well-known in the art. Chips of early code 69B (FIG. 2D) start and stop O.SAtemi 
seconds before the corresponding chips of prompt PRN code 29 (FIG. 2D), and the 
chips of late code 69C (FIG, 2D) start and stop O.SAtemi seconds after prompt PRN 
code 29 (FIG, 2D). Early-minus-late PRN code 35 (FIG. 2D) is the difference 
between early code 69B (FIG. 2D) and late code 69C (FIG. 2D). Example segments 
of these four types of replica codes are depicted in FIG. 2D. 
[00060] Referring to FIGs. 1, 2A, 2C, and 2D software correlator 19 (FIG. 1) 
receives, through conventional RF front end 13 and DAQ system 17, signal data 21 
(FIG. 1), the raw data 12 (FIG. 1) source of which is sampled at the rate fs = \IAts Hz. 
In order to process the resulting RF signal data 21, software correlator 19 (FIG. 1) 
needs prompt PRN code 29 (FIG. 2D) and early-minus-late PRN code 35 (FIG. 2D) 
replicas sampled at the same times as raw signal 12 (FIG. 1). FIG. 2D depicts si?cteen 
sample times as vertical dash-dotted lines. Referring to FIG. 2D, prompt PRN code 29 
(FIG. 2D) can be represented by its prompt PRN code sign 29A (FIG. 2D) at the 
sample times. The bit value one represents +1, and the bit value zero represents -1 . 
Prompt PRN code sign 29A (FIG. 2D), shown at the sixteen sample times — starting 
with three Is, continuing with ten Os, and finishing with another three Is — is a 16-bit 
word stored as the integer 2^^+2^V2^^+2V2^+2^ = 57351. Early-minus-late PRN code 
35 (FIG. 2D) requires a 1 .5-bit representation. A zero mask bit is set to zero if early- 
minus-late PRN code 35 takes on the value zero, and it is set to one if early-minus-late 
PRN code 35 equals +2 or -2, Early-minus-late PRN code zero mask 35B (FIG. 2D) 
at sixteen sample times shown in FIG. 2D is equivalent to 2^^+2^V2^ = 12292. A 2's 
sign bit is set to one if early-minus-late PRN code 35 (FIG. 2D) equals +2 at the 
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sample time, and it is set to zero if the code equals -2, The 2's sign bit is irrelevant if 
the corresponding early-minus-late PRN code zero mask 35B (FIG. 2D) bit equals 
zero. Early-minus-late PRN code sign 35A (FIG. 2D) for sixteen sample times 
contains X values that indicate bits whose values are irrelevant because the 
corresponding early-minus-late PRN code zero mask 35B (FIG. 2D) bits are zero. In 
an illustrative embodiment, all the X values become zero, thus the equivalent integer 
for early-minus-late PRN code sign 35A (FIG. 2D) is 2^ = 4. 

[00061] Continuing to refer to FIG. 3 A, an alternative to taking the prompt PRN 
code 29 and early-minus-late PRN code 35 from pre-computed over-sampled PRN 
code table 28 is to generate prompt PRN code sign 29 A, early-minus-late PRN code 
sign 3 5 A, and early-minus-late PRN code zero mask 35B using real-time over- 
sampled PRN code generator 30A (FIG. 3A). Shown in FIG. 3 A are two circles and a 
loose arrow with a quarter circle pointer. These are the symbols for a switch and 
indicate the ability of the system to choose possible alternate sources of PRN code. 
Using the real-time over-sampled PRN code generator 30A includes a step of 
generating the PRN code chips in real-time by conventional means. For example, the 
GPS civilian L2 CL and CM codes are generated by a 27-bit feedback shift register 
(see The New L2 Civil Signal, R.D. Fontana et al.. Proceedings of the ION GPS 200 L 
September 1 1-14, 2001, Salt Lake City, UT, pp. 617-631). The method further 
includes the steps of choosing chip values from the PRN code, where the chip values 
correspond to a data interval that contains the samples of a data word and where the 
chips have a known timing relative to the data interval, transforming the relative 
timing into a time grid index, and translating the PRN code chip values and the time 
grid index for the data interval into the PRN code's over-sampled bit-wise parallel 
format. These latter steps can be carried out efficiently by using a table look up 
function. One table each for prompt PRN code sign 29 A, early-minus-late PRN code 
sign 3 5 A, and early-minus-late PRN code zero mask 35B can include integer values 
that constitute the bit-wise parallel representation of the PRN code for the sample 
times associated with the data word in question. Indices into each 1 -dimensional table 
are functions of (a) the time offset between the first PRN code chip and the first 
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sample time of the given data word, and (b) the bit pattern of the PRN code chips that 
span the sample times of the data word. The sizes of the tables are independent of the 
period of the PRN code that is being over-sampled. The tables can be re-used for 
multiple PRN codes in a multi-channel receiver. The computation and use of the 
tables are discussed in more detail later. 

[00062] Continuing with the description of bit-wise parallelism with respect to 

the operations of software correlator 19 (FIG. 1), and continuing to refer to FIG. 3A, 
the specially-developed algorithms described herein make use of bit-wise parallelism 
so that a single programming language statement, such as a C code command, can 
partially process up to thirty-two samples at a time. Previously referred-to carrier 
replica signal 25 in the form of cosine and sine signals are stored as binary carrier 
replica sign 25A and carrier replica magnitude 25B. The format of this representation 
is defined in Table 2 and illustrated in FIG. 3C which is a reconstructed carrier and 
carrier replica in the form of representative sine signal 51 (FIG. 3C) shown in optimal 
2-bit representation 53 (FIG. 3C) that has the minimum square error. The format of 
Table 2 assumes that the cosine and sine signals have an amplitude of approximately 
2.4. Note that other representations, beyond 2-bit representation, are possible. In 
general, more bits yield a better SNR, but can also require a larger number of 
computations for the correlation operations. 



Carrier Replica Sign 
25A 


Carrier Eeplica 
Magnitude 2SB 


Carrier Replica Value 


0 


0 


-1 


0 


1 


-2 


1 


0 


+1 


1 


1 


+2 



Table 2. Sign and magnitude combinations of the stored intermediate-frequency 
baseband mixing carrier sine wave replicas and the values that they represent 



[00063] Continuing to refer to FIG. 3 A, multiplication of the RF front-end 
output representation, the signal sign 21 A and signal magnitude 2 IB, of Table 1 by the 
sine wave representation, carrier replica sign 25A and carrier replica magnitude 25B, 
of Table 2 yields baseband mixed signals 23, consisting of baseband mixed sign 23 A, 
baseband mixed high magnitude 23B, and baseband mixed low magnitude 23 C, that 
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can take on the values -6, -3, -2, -1, +1, +2, +3, and +6, as shown in Table 3. 
Baseband mixed high magnitude 23B is simply signal magnitude 2 IB, and baseband 
mixed low magnitude 23C is carrier replica magnitude 25B. Thus, these two 
magnitude bits are available without the need for computation. Baseband mixed sign 
23 A is the result of an EXCLUSIVE OR operation between signal sign 21 A and 



carrier replica sign 25 A. Notice how the relationship of the sign bit value with the 
actual sign gets reversed from that of Tables 1 and 2. 



Baseband 
Mixed Sign 
23A 


Baseband Mixed 
High Magnitude 
23B 


Baseband Mixed 
Low Magnitude 
23C 


Baseband Mixed 
Value 


0 


0 


0 


+1 


0 


0 


1 


+2 


0 


1 


0 


+3 


0 


1 


1 


+6 


1 


0 


0 


-1 


1 


0 


1 


-2 


1 


1 


0 


-3 


1 


1 


1 


-6 



Table 3. Sign, high-magnitude, and low-magnitude combinations of the baseband- 
mixed signal and their corresponding signal values. 



[00064] Continuing to refer to FIG. 3 A, and continuing to describe the bit-wise 
parallel algorithms, the required amount of storage for tables of pre-computed prompt 
PRN code 29 and early-minus-late PRN code 35 can be greatly reduced by making 
two simplifications. First, the prompt PRN code 29 is stored as prompt PRN code sign 
29 A. This representation is shown in Table 4. The early-minus-late PRN code 35, on 
the other hand, is stored in a two-bit representation (actually a 1.5 bit representation): 
early-minus-late PRN code sign 35A and early-minus-late PRN code zero mask 35B, 
as denoted in Table 5. Note that the X in the first column of Table 5 indicates that 
zero or one can be placed in this location without affecting the corresponding code 
value. The X signifies a lack of effect of the sign bit on the code value when the zero 
mask bit equals zero. This is why the early-minus-late PRN code 35 representation is 
referred to as a 1.5-bit representation. This X value will aflFect the corresponding fiilly 
mixed early-minus-late integrand sign 33 A, but it will not affect any of the early- 
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minus-late value words because the zero value in the corresponding zero mask 
location will null out the corresponding bit of all early-minus-late value words. 



Promptt PRN Code 
Sign 29A 


Prompt Code Value 


1 


+1 


0 


-1 



Table 4. Sign bits of the prompt C/A code and the corresponding prompt signal values. 



Early-imiDiniiuis-late PRN 


Early-minus-Date PMN 


Eariy-Minns-Laite Code 


Code Sign 3SA 


Code Zero Mask 3SB 


Value 


X 


0 


0 


0 


1 


-2 


1 


1 


+2 



Table 5. Sign and zero mask bit combinations of the early-minus-late PRN code 35 
and the corresponding signal values. 



[00065] Another simplification in the pre-computed over-sampled PRN code 
table 28, and continuing to refer to FIG. 3 A, can be to ignore code Doppler shift 
variations. All signals in the table are assumed to have zero Doppler shift; i.e., all C/A 
codes in the table assume that fyt+i -Tjk = 0.001 sec. Note that the period of 0.001 is 
applicable for accumulations that use the fall 1023 chips of the C/A code only. Any 
other type of code or accumulation interval may have a different period. The code 
phase errors due to this assumption can be eliminated by choosing a replica code fi-om 
the pre-computed over-sampled PRN code table 28 whose midpoint occurs at the 
desired midpoint time (rjk + TjmV^ . The only other effect of this assumption can be a 

small correlation power loss, which is no more than 0.014 dB if the magnitude of the 
Doppler shift is less than 10 KHz. The pre-computed over-sampled PRN code table 
28 should include a selection of different phases, for example fourteen, as measured 
relative to a signal sample spacing of, for example, 175 nsec. This translates into a 
code phase spacing of, for example, 12.5 nsec, which equals a pseudo range 
measurement digitization level of 3.8m, or a maximum measurement error of 1.9m. 
The number of phases in the pre-computed over-sampled PRN code table 28 is 
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dependent upon the design of the system and no set number of phases is required by 
the present invention. Referring to FIG 6, suppose that pre-cx)mputed over-sampled 
PRN code table 28 stores over-sampled bit-wise parallel representations of chips C(l) 
through C(A4), The table must allow for the retrieval of over-sampled bit-wise parallel 
code replicas for a range of start times of code chip C(l) that span the entire first data 
sample word in the accumulation interval Wi 95 (FIG 6). The table may contain code 
replicas whose different phases )deld start times that span only a single sample interval 
of data word W\ 95 (FIG 6), which is only l/w, of the required number of start times. 
In this case the software correlator may apply bit shift operations to a tabulated PRN 
code replica fi'om that sample interval in order to generate the over-sampled bit-wise 
parallel PRN code replica that applies when chip C(l) starts in a different sample 
interval of data word Wi 95 (FIG 6). 

[00066] Continuing to refer to FIG. 3 A, and further continuing to describe the 
bit- wise parallel algorithms, prompt PRN code 29 and early-minus-late PRN code 35 
replicas can be mixed with the baseband mixed signals 23 to form fully mixed prompt 
integrand 3 1 by an EXCLUSIVE OR operation and bit re-definitions. An 
EXCLUSIVE OR between prompt PRN code sign 29A and baseband mixed sign 23 A 
produces fiilly mixed prompt integrand sign 31 A given in Table 6. The fiilly mixed 
prompt integrand high magnitude 3 IB and fiilly mixed prompt integrand low 
magnitude 3 IC are baseband mixed high magnitude 23B and baseband mixed low 
magnitude 23C, also given in Table 6. Note that the Table 6 representation is identical 
to that of Table 3 except for the inversion in the meaning of the sign bits. The number 
of magnitude bits is dependent upon the design of the system and no set number of 
magnitude bits is required by the present invention. A change in the number of 
magnitude bits will cause a change in the number of entries of the equivalent of Table 
6 and it will affect the possible values of the integrand. 
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FuUv Mixed 
PromDt 
Integrand Sign 
31A 


Fiillv Mix^d Prnmnt 
Integrand Hi^h 
Magnitude 31B 


Fully Mixed 
PromDt 
Integrand L w 
Magnitude 31C 


FuUv Mixed 
PromDt 
Integrand 
Value 


0 


0 


0 


-1 


0 


0 


1 


-2 


0 


1 


0 


-3 


0 


1 


1 


-6 


1 


0 


0 


+1 


1 


0 


I 


+2 


1 


1 


0 


+3 


1 


1 


1 


+6 



Table 6. Sign, high-magnitude, and low-magnitude bit combinations of the fully 
mixed prompt integrand 31 and its corresponding values. 



[00067] Still continuing to refer to FIG. 3 A, the mixing of the early-minus-late 
PRN code 35 with the baseband mixed signals 23 forms ftilly mixed early-minus-late 
integrands 33. Fully mixed early-minus-late integrand sign 33A is an EXCLUSIVE 
OR between early-minus-late PRN code sign 35A and baseband mixed sign 23 A. 
Fully mixed early-minus-late integrand high magnitude 33B and fiiUy mixed early- 
minus-late integrand low magnitude 33C are, as above, baseband mixed high 
magnitude 23B and baseband mixed low magnitude 23C. Fully mixed early-minus- 
late integrand zero mask 33D is early-minus-late PRN code zero mask 35B. The 
resulting representation is given in Table 7. As in Table 5, each X entry in the table 
indicates that the corresponding bit can be either zero or one without affecting the 
corresponding integrand value. 
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Table 7. Sign, high-magnitude, low-magnitude, and zero mask bit combinations of the 
fiilly mixed early-minus-late integrands 33 and their corresponding values. 



[00068] Referring now to FIGs. 3 A, 3B, 4A, and 4B, the method for computing 
in-phase and quadrature accumulations for every accumulation period, for example 
every millisecond for GPS C/A code, by use of bit-wise parallelism includes the steps 
of selecting carrier replica signal 25 (FIG. 3 A) according to the proximity of its 
frequency to the desired frequency, and representing sample signal data 21 (FIG. 3 A) 
and carrier replica signal 25 (FIG. 3 A) from at least one channel as bits in signal sign 
21 A (FIG. 3A) and, if present, signal magnitude 21B (FIG. 3A) and carrier replica 
sign 25 A (FIG. 3 A) and carrier replica magnitude 25B (FIG. 3 A) (method step 101, 
FIG. 4A). Note that carrier replica signal 25 (FIG. 3 A) is chosen so that its frequency 
is close to the correct signal frequency. The method also includes the step of mixing 
signal data 21 (FIG. 3 A) to baseband by computing in-phase and quadrature baseband 
mixed sign 23 A (FIG. 3 A) and in-phase and quadrature baseband mixed high and low 
magnitude 23B/C (FIG. 3A) (method step 103, FIG. 4A), The method further 
includes the steps of selecting PRN code from pre-computed over-sampled PRN code 
table 28 (FIG. 3 A) or of computing it using real-time over-sampled PRN code 
generator 30A (FIG. 3 A), representing prompt PRN code 29 (FIG. 3 A) as prompt 
PRN code sign 29 A (FIG. 3 A), and representing early-minus-late PRN code 35 (FIG. 
3A) from as early-minus-late PRN code sign 35A (FIG. 3A) and early-minus-late PRN 
code zero mask 35B (FIG. 3 A) (method step 105, FIG. 4 A). The method further 
includes the step of de-spreading in-phase and quadrature baseband mixed signal 23 
(FIG. 3 A) by mixing it with prompt PRN code 29 (FIG, 3 A) and early-minus-late 
PRN code 35 (FIG. 3 A), resulting in in-phase and quadrature fiiUy mixed prompt 
integrands 31 (FIG. 3 A), and fully mixed early-minus-late integrands 33 (FIG. 3 A) 
(method step 107, FIG. 4A). The method further includes the step of using prompt 
value word logic 27A (FIG. 3B) to compute prompt integrand value words 27 (FIG. 
3B) from the in-phase and quadrature fiiUy mixed prompt integrands 3 1 (FIG. 3 A). 
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The method fiirther includes the step of using early-minus-late value word logic 37 A 
(FIG. 3B) to compute early-minus-late integrand value words 37 (FIG. 3B) from the 
folly mixed early-minus-late integrands 33 (FIG. 3B) (method step 109, FIG, 4A). 
The method further includes the steps of summing over each prompt integrand value 
word 27 and early-minus-late integrand value word 37 (FIG. 3B) the number of one 
bits (or zero bits) using one bits summation table 38 (FIG. 3B) or using a processor 
command if available (method step 111, FIG. 4B), and summing, over the 
accumulation interval, the number of one bits (or zero bits) in each prompt integrand 
value word 27 and early-minus-late integrand value word 37 to produce prompt 
accumulations 41 (FIG. 3B) and early-minus-late accumulations 49 (FIG. 3B) (method 
step 113, FIG 4B). The method forther includes the step of multiplying prompt 
accumulations 41 (FIG. 3B) and early-minus-late accumulations 49 (FIG. 3B) by 
corresponding values 41 A and summing the results over the value words of each 
signal for an entire accumulation interval to yield in-phase and quadrature summed 
prompt accumulations 45 (FIG. 3B) and sununed early-minus-late accumulations 47 
(FIG. 3B) (method step 1 1 5, FIG. 4B) that are stored for use by acquisition techniques 
or tracking loops. The method forther includes the step of rotating the in-phase and 
quadrature summed prompt accumulations 45 (FIG. 3B) and summed early-minus-late 
accumulations 47 (FIG. 3B) (method step 117, FIG. 4B) to simulate a condition in 
which baseband mixing had been performed using cosine and sine signal replicas with 
the correct frequency and phase. If there are more channels to process (decision step 
119, FIG. 4B), the method includes the step of repeating the previous steps beginning 
at method step 101, FIG. 4A. If there are no more channels to process (decision step 
119, FIG. 4B), the method includes the step of setting parameters for the next 
accumulation period, including storing current C/A code phases, epoch counters, 
carrier phases, and carrier Doppler shifts (method step 121, FIG. 4B). If the time 
period to wait until the next accumulations need to be calculated has not expired 
(decision step 123, FIG. 4BX the method includes the step of sleeping until the 
expiration of the time period (method step 125, FIG. 4B). If the time period has 
expired (decision step 123, FIG. 4B), the method includes the step of repeating the 
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previous steps beginning at method step 101, FIG. 4 A, The length of the time period 
depends on the nominal accumulation period. It is set to be less than this period, 
normally between 50% to 90% of this period, to reduce the possibility that 
accumulations are missed for any channels. 

[00069] Referring again to FIGs. 3B and 4A, method step 109 (FIG. 4A) calls 
for computing value words. This computation starts by performing bit-wise parallel 
Boolean logic for each of the possible values in the right-hand column of the prompt 
integrand representation in Table 6. A 32-bit prompt integrand value word 27 (FIG. 
3B) is computed for each thirty-two samples and each row of Table 6, The prompt 
integrand value word 27 (FIG. 3B) contains ones for the sample times when the actual 
integrand equals the corresponding value in the right-hand column of Table 6, and 
zeros for the remaining times when the actual integrand does not equal this value. The 
prompt integrand value words 27 (FIG. 3B) corresponding to the possible Table 6 
values are formed by method step 109 (FIG. 4A) as follows: 

MBSfUSONE = NOT(SIGN) AND [NOT(HIGHMAG) AND NOT(LOWMAG)] (13) 



MDSrUSTWO = NOT(SIGN) AND [NOT(HIGHMAG) AND LOWMAG] (14) 

MNUSTHREE = NOT(SIGN) AND [fflGHMAG AND NOT(LOWMAG)] (15) 

MINUSSIX = NOT(SIGN) AND [fflGHMAG AND LOWMAG] (16) 

PLUSONE = SIGN AND [NOT(fflGHMAG) AND NOT(LOWMAG)] (17) 

PLUSTWO = SIGN AND [NOT(fflGHMAG) AND LOWMAG] (18) 

PLUSTHREE = SIGN AND [fflGHMAG AND NOT(LOWMAG)] (19) 

PLUSSIX = SIGN AND [fflGHMAG AND LOWMAG] (20) 



[00070] Continuing to refer to FIGs. 3 A, 3B, 4 A, and 4B, method steps 109 
(FIG 4 A), 1 1 1 (FIG. 4BX and 1 13 (FIG. 4B) call for operations for the folly mixed 
early-minus-late integrands 33 (FIG. 3A) that are similar to those for the folly mixed 
prompt integrands 31 (FIG. 3 A). Early-minus-late integrand value words 37 (FIG. 
3B) correspond to values that are double those of the prompt integrand value words 27 
(FIG. 3B), i.e., the MINUSSIX word becomes the MINUSTWELVE word. Also, an 
additional AND operation must be performed with the zero mask bits of Table 7 in 
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order to mask out sample times when the early and late PRN codes cancel each other. 
Possible formulas for the method step 109 (FIG. 4A) computation of these early- 
minus-late integrand value words 37 (FIG. 3B) are as follows: 
MINUSTWO = [ZEROMASK AND NOT(SIGN)] AND [NOT(HIGHMAG) AND 

NOT(LOWMAG)] (21) 
MINUSFOUR = [ZEROMASK AND NOT(SIGN)] AND [NOT(fflGHMAG) AND 

LOWMAG] (22) 
MINUSSIX = [ZEROMASK AND NOT(SIGN)] AND [fflGHMAG AND 

NOT(LOWMAG)] (23) 
MDSfUSTWELVE = [ZEROMASK AND NOT(SIGN)] AND [fflGHMAG AND 

LOWMAG] (24) 
PLUSTWO = [ZEROMASK AND SIGN] AND [NOT(fflGHMAG) AND 

NOT(LOWMAG)] (25) 
PLUSFOUR = [ZEROMASK AND SIGN] AND [NOT(fflGHMAG) AND 

LOWMAG] (26) 
PLUSSIX = [ZEROMASK AND SIGN] AND [fflGHMAG AND 

NOT(LOWMAG)] (27) 
PLUSTWELVE = [ZEROMASK AND SIGN] AND [HIGHMAG AND 

LOWMAG] (28) 
Additional zero masking can occur in the first and last words of an accumulation 
interval. This is true because the start and stop times of an accumulation interval do 
not normally fall at the boundaries of data words. Therefore, the bits in the first word 
that precede the accumulation interval may need to get zero masked as might the bits 
in the last word that come after the end of the accumulation interval. 
[00071] Referring primarily to FIGs. 3B and 4B, the one bits counting 
operations of method step 111 (FIG. 4B) form the count of the number of one bits in 
each of the eight value words. If there are no such counting operations in the 
instruction set of microprocessor 16 (FIG. 1), the counting can be accomplished using 
a table look-up. In the case of a table look-up, prompt integrand value words 27 and 
early^minus-late integrand value words 37 (FIG. 3B) can be used as addresses in one 
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bits summation table 38 (FIG. 3B), and one bits summation table 38 (FIG. 3B) can 
output the number of one values (or zeros) in the address. For example, if the table 
look-up operation is called BITSUM, the following computations can be performed to 

compute one-bits counts: 

ONESCOUNT = BITSUM(VALUEWORD) (29) 
where the output of the table ONESCOUNT is the number of one bits in the word 
VALUEWORD, This operation is repeat^ for each of the prompt integrand value 
words 27 (FIG. 3B) and early-minus-late integrand value words 37 (FIG. 3B) in order 
to accomplish method step 111 (FIG. 4B). Selection of table width, for example 16- 
bit or 32-bit, depends on the amount of memory available and other design decisions. 
If the table width is smaller than the number of bits in a value word, then multiple 
calls of the table are used in order to sum up the total number of one values in a given 
value word. Each call takes as input only a portion of the bits in the value word. 
[00072] Continuing to refer primarily to FIGs. 3B and 4B, the accumulation 
operations of method steps 1 13 (FIG. 4B) and 1 1 5 (FIG. 4B) sum the one bit counts 
for each prompt integrand value word 27 (FIG. 3B) and for each early-minus-late 
integrand value word 37 (FIG. 3B) over the entire accumulation interval, multiply 
each result by the value 41 A (FIG. 3B) that is associated with the value word, and 
sum all of these scaled value accumulations to form the accumulations of equations (6) 
and (7), summed prompt accumulation 45 (FIG. 3B) and sunmied early-minus-late 
accumulation 47 (FIG. 3B). For example, the following computations can be 
performed to compute the in-phase summed prompt accumulation 45 in equation (6) 
as follows: 



Ig/ki^) = - E ONESCOUNT(-l)/^/ -2 Z ONESCOUNT(-2)/^/ 



-3 Z 0NESC0UNT(-3)/^/ - 6 Z 0NESC0UNT(-6)/ / 
/=1 /=1 
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+ Z ONESCOUNT(+l)/p/ + 2 Z 0NESC0UNT(+2)/./ 
/=! /=1 

+ 3 2: 0NESC0UNT(+3) /«/ + 6 Z ONESCOUNT(+6) /„/ (30) 

/=l /=! 

where / is the index of successive bit-wise parallel data words in the accumulation 
interval, is the total number of data words in the interval, and ONESCOUNT(A:)/p/ 
is the ones count for the corresponding value word 41 (FIG. 3B) associated with value 
A: 41 A (FIG. 3B) for the Ah data word interval and the in-phase summed prompt 
accumulation 45 (FIG. 3B). The quadrature summed prompt accumulations 45 (FIG. 
3B) and the in-phase and quadrature summed early-minus-late accumulations 47 (FIG. 
3B) are calculated in a similar manner. The only difference is in the actual 
ONESCOUNT values used and, for the case of early-minus-late signals, the set of k 
values 41 A (FIG. 3B). 

[00073] Continuing to refer primarily to FIGs. 4A and 4B, the method of the 
present invention can be adapted to work with a different number of bits in the 
representation of the RF front-end output and of the baseband mixed signals. An 
increase above two bits can make the logic more complex and may decrease the time 
savings over str^ght integer arithmetic. A decrease to a 1-bit representation can have 
the opposite effect. For example, if the RF front-end uses 1-bit digitization rather than 
2-bit digitization while carrier replica signal 25 (FIG. 2B) retains its 2-bit digitization, 
then the operation count can decrease by a factor of almost two for the 1-bit method, 
which can make the logic execute about 4.2 times faster than straight integer 
arithmetic. 

[00074] Returning to the discussion of determining PRN code, and now 
referring again FIGs. 2C, 2D, and 3 A, the real-time generation of bit-wise parallel 
over-sampled prompt PRN code sign 29A (FIGs. 2D and 3 A), early-minus-late PRN 
code sign 35A (FIGs. 2D and 3A), and early-minus-late PRN code zero mask 35B 
(FIGs. 2D and 3 A) can be carried out by real-time over-sampled PRN code generator 
3 OA (FIG. 3 A). The inputs to this calculation are the actual PRN code chip length 65 
(FIG. 2D), Ate, the sample interval 63 (FIG. 2D), Als, the nominal early-to-late code 
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delay 61 (FIG. 2D), Atemh the end time of the first code chip relative to the first sample 
time, or put another way, the time lag Aio, 67 (FIG. 2D) fi^om the first RF sample time 
to the end time of the first prompt PRN code chip, and prompt code chips 91 (FIGs. 
2D and 6). The outputs are the three integers that store the prompt PRN code sign 
29A (FIGs. 2D and 3 A), early-minus-late PRN code zero mask 35B (FIGs. 2D and 
3A), and early-minus-late PRN code sign 35A (FIGs. 2D and 3A), which are all in bit- 
wise parallel format. 

[00075] Referring again to FIGrs. 2C and 2D, table look-ups can be used to 
translate a PRN code and its timing information to bit-wise parallel representations of 
its over-sampled prompt and early-minus-late versions. The required table look-ups 
can be simplified by recognizing that the following parameters are substantially 
constant, for the purposes of this calculation: sampling interval 63 (FIG. 2D), Ats, the 
nominal chip length , Atcnom, the early-minus-late code delay 61 (FIG. 2D), Atemu used 
by software correlator 19 (FIG 1), and the maximum number of chips that span a data 
word of microprocessor 16 (FIG. 1). The difference between the actual chipping rate 
4/c (reciprocal of Ate) and the nominal chipping rate ^cnom (reciprocal of Atcnom) that is 
used for the above simplification can be accommodated by correcting time lag 67 
(FIG. 2D), Ato, for the average effects of Doppler shift, a procedure discussed later. 
Using the simplification, each look-up table has two variable inputs: the actual set of 
prompt code chips 91 (FIGs. 2D and 6) and time lag 67 (FIG. 2D), Ato. A table look- 
up procedure for each signal component yields a single integer result for prompt PRN 
code sign 29A (FIGs. 2D and 3 A), another single integer result for early-minus-late 
PRN code zero mask 35B (FIGs. 2D and 3 A), and yet another single integer result for 
early-minus-late PRN code sign 35 A (FIGs. 2D and 3 A). 

[00076] Time lag 67 (FIG. 2D), Ato, can take on an infinite number of values in 
the continuous range: 

-\At,^i <Ato<At, -\At,^i (3 1) 

This range's lower limit guarantees that the end time of the first late chip occurs no 
earlier than the first sample time. A lower time lag 67 (FIG. 2D) Alo value would 
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make the first chip irrelevant to the prompt PRN code 29 (FIG. 2D), early code 69B 
(FIG. 2D), and late code 69C (FIG. 2D) at all of the sample times. The upper limit in 
equation (31) guarantees that the start time of the first late chip occurs no later than the 
first sample. A larger value of Ato would leave the late code 69C (FIG. 2D) at the first 
sample time undefined based on the available code chips. 

[00077] Referring now to FIG. 5, to create an electronically processable table, 
the continuous range of ^d/o values can be replaced with a discrete grid having m 
equally spaced points per sample interval 63 (FIG. 2D), Ats. The integer m is chosen 
to be large enough so that the granularity Ats/m gives sufficient PRN code timing 
resolution. In GPS applications m is usually chosen to be large enough so that 
(cAtJm) is on the order of several meters or less, where c is the speed of light, but 
reasonably sized because the table sizes are usually proportional to m. Given a choice 
of m, the grid of relative end times of the first prompt code period is: 
kAt 

Alo, = — ^ ioTk = k„^,...,k„^ (32) 
m 

where the limits 

*,,„=yZoor(-^^^).2 (33a) 

kmax=floor{ ) (33b) 

provide full coverage of the interval defined in equation (31). The floor() function 
rounds to the nearest integer in the direction of -oo. This kmin value can cause the 
minimum A/^j^to fall slightly below the lower limit in equation (3 1), which can cause 

memory inefficiency, but this value is advantageous because it may simplify some 
further computations. 

[00078] The size for each table can be a function of the maximum number of 
code chips that may fall within a data word's sample range. Given Ato^ bit information 
for the following number of code chips is required in order for the prompt PRN code 
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29 (FIG. 2D), early code 69B (FIG. 2D), and late code 69C (FIG. 2D) to be ftilly 
specified at all of the data word's sample times: 



liAto) = floor 



► + 2 (34) 



where rts is the number of data samples that can be stored in bit-wise parallel format in 
each word. It is clear from equation (34) that /(zl/o) is a non-increasing function of Mq. 
Therefore, the maximum number of required chips occurs at the minimum value of 

AtQ. 

^ = /(^o*»*J (35) 
The size of each table can be determined from the parameters kmm, kmax, and X. The 
grid contains ktot = (^max-^m/w+l) different time offsets of the first code chip. At each of 
these grid points there are 2^ possible combinations of the code chips. Thus, each 
table optimally contains ktotxT^ entries, and each entry is optimally an unsigned integer 
in the range from 0 to 2"'-l . 

[00079] Continuing to refer to FIG. 5, each table can be stored as an array with 
a single index. The first 2^ entries correspond to the 2^ different possible chip 
sequences that can occur at Ato = ^oAnmn > entries correspond to Ato = 

^0(A7w/>i+l) > so forth. The tabulated bit sequences for a fixed Ato are ordered by 

interpreting the sequence as a binary index counter with the first chip being the most 
significant counter bit and the chip being the least significant bit. The integer 
elements of the table can be the x{i) table elements 81 with corresponding code time 
offset 83 Atok, and corresponding bit sequence 85 of the chips. The array index of a 
given x(/) table element 81 can be computed based on its code time offset 83 Atokgrid 
index k and its corresponding bit sequence 85. The corresponding bit sequence 85 
associated with the array index consists of the chip values C(l), C(2), C(3), ... C(L). 
The C{j) chip values are either zero or one, with zero representing a -1 PRN code 
value and one representing a + 1 PRN code value, and they are listed in order of 
increasing time. The corresponding array index of the x{i) table element 81 is: 
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/[^.C(1),C(2).C(3),...,C(L)] = \ + (k- k^i„) X 2^ + l.CU)2^~j 

7=1 

This equation can be inverted to give the code time offset 83 grid index k and the 
corresponding bit sequence 85 as functions of the x(i) table element 81 index /: 

kii) = k^i„+floor[(i-l)/2^] 



(36) 



CO",/) = mod 



floor 



fmod[0 -1),2'^]1 



1 



,2 



for>=l,2,3,...,£ 



(37a) 
(37b) 



where mod(y,z) =y- zxfloor(y/z) is the usual remainder function. 
[00080] Continuing to refer to FIG. 5, the following computations generate the 
x(i) table elements 81 entries of the three tables. Given i, the corresponding code time 
offset 83 grid index k(i) is computed from equation (37a) and is used to generate three 
sequences of chip indices: 




n-l- 



k(i) 




1 




m 








k(i)^ 

m 




+ 




k(i)-] 








m \ 






2At, 



forw= 1,2, 3, ...,ns 
for«= 1,2, 3, «5 
for«= 1, 2, 3, «i 



(38a) 
(38b) 
(38c) 



where n is the index of the sample time within the over-sampled data word. The 
integer jp(n,i) is the index of the PRN code chip that applies at sample n for the prompt 
PRN code 29 (FIG. 2D). The integers jeO^J) and j^nj) are defined similarly for the 
early code 69B (FIG. 2D) and late code 69C (FIG. 2D), respectively. The formulas in 
equations (38a)-(38c) amount to time measurements of each sample given in units of 
chip lengths past the first chip. These indices, in turn, can be used to determine the 
chip values that apply at the sample times: 

C,{nJ) = C[j,(nJ);i] for /i = 1, 2, 3, (39a) 
Q(«.0 = CUMi);i] for /I = 1, 2, 3, (39b) 
^(''^O = C[j, (nj);i] for « = 1, 2, 3, (39c) 
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where Cp{n,i) is the over-sampled prompt PRN code 29 (FIG. 2D), and Ce{n,i) and 
Cin,i) are, respectively, the early code 69B (FIG. 2D) and late code 69C (FIG. 2D). 
Each of these code bit values is either zero or one, as dictated by the outer mod{ ,2) 
operation in equation (37b). These over- sampled chip values can, in turn, be used to 
formulate tabulated functions Xp(/), Xemizmij), and Xemas{i) that generate the unsigned 
integers that constitute the bit-wise parallel code representations of the three tables: 

^p(0 = |;C,(«,/)x2'''-" (40a) 
Xemlzmii) = Zmod{{C,inJ) + C/(/i,i)].2} X 2'^-"" (40b) 

^emllsii) = Zmod{[C,{nJ) + C/(/i,i)],2} x C,{nJ) x 2"-"'' (40c) 

where Xp(/) is the entry of the prompt sign table, Xemizmif) is the entry of the early- 
minus-late zero mask table, and Xemnsif) is the entry of the early-minus-late 2*s sign 
table. Note that the formula used in equation (40c) is only an example illustrative 
embodiment of the early-minus-late 2*s sign table calculation. It places zeros in all of 
the X entries of early-minus-late PRN code sign 35A (FIG. 2D). There exist alternate 
formulas that are equally correct but that do not place zeros in the X entries. 
[00081] The table layout in FIG. 5 is only an illustrative embodiment of how 
one can construct a table that can be used to translate PRN code chip values and 
timing information into data words that store the bit-wise parallel representations of 
the over-sampled prompt PRN code sign 29A (FIG. 2D), early-minus-late PRN code 
zero mask 35B (FIG. 2D), and early-minus-late PRN code sign 35A (FIG. 2D). Other 
table layouts are also possible. Possible illustrative index calculations are described 
below for indexing into the tables for PRN code retrieval during accumulation 
calculations If another table layout is used, then different indexing calculations might 
be needed. Furthermore, different indexing calculations can be used even for the 
illustrative table layout shown in FIG. 5. 

[000S2] Referring now primarily to FIG. 6, accumulation calculations, as have 
been previously outlined herein and elsewhere, work with a fixed sequence of code 
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chips. The prompt version of this sequence has a specified timing relationship to the 
incoming RF signal data 21 (FIG. 2 A). This relationship can be pre-determined by a 
code search algorithm if software receiver 10 (FIG. 1) is in acquisition mode or by its 
delay-locked loop if it is in tracking mode. Software correlator 19 (FIG. 1) can 
calculate an accumulation using prompt code chips 91 (FIG. 6) C(l) through C(M), 
The timing of the prompt replicas of prompt code chips 91 (FIG. 6) can define the 
accumulation interval. The chip sequence starts at start lag 93 (FIG. 6) jdistart seconds 
past the first sample of data word W\ 95 (FIG. 6), it chips at the constant chipping rate 
fc = l/Atc, and it ends at end time 97 (FIG. 6), which occurs Atstar&MAtc seconds after 
the first sample of data word Wi 95 (FIG. 6). The end of the A/* prompt code chip can 
occur during data word Wn 99 (FIG. 6), which implies that 

N = ceill ^ (41) 

I ^s^s J 

where the ceilQ fiinction rounds to the nearest integer towards -i-oo. Some of the initial 
bits of data word Wi 95 (FIG. 6) and some of the final bits of data word Wm 99 (FIG. 
6) may not be included in the accumulation. Let riexo be the number of initial bits of 
data word Wi 95 (FIG. 6) that are excluded, and let rtexfht the number of final bits of 
data word Wn 99 (FIG. 6) that are excluded. The timing relationship in FIG. 6 
implies that these numbers are: 



= ceil 



V J 



n^^ ^ nN - ceil 



Ai, 



(42a) 
(42b) 



[00083] These sample counts can be used to develop additional zero mask 
words that software correlator 19 (FIG. 1) uses to properly process the first and last 
data words during its bit-wise parallel accumulation calculations, as defined in^ 72- 
Channel Real-Time GPS LJ Software Receiver, B. M. Ledvina et al.. Proceedings of 
the ION National Technical Meeting . January 22-24, 2003, Anaheim, CA. and Bit- 
Wise Parallel Algorithms for Efficient Software Correlation Applied to a GPS 
Software Receiver. B.M. Ledvina et al., to appear in the IEEE Transactions on 
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Wireless Communications, 2003, both incorporated herein in their entirety by 
reference. Note that equations (41)-(42b) and all related timing considerations herein 
use the following code chip start/stop convention: a sample is correlated with a 
particular code chip if the start time of the code chip coincides exactly with the sample 
time, but it will not get correlated with that chip if its sample time coincides exactly 
with the end time of the code chip. 

[00084] Continuing to refer to FIG. 6, efficiently determining the correct Xp(i\ 
Xemtm(iX and Xemns(i) bit-wise parallel code representations for the ^ data words Wi 95 
(FIG. 6) through data word Wn 99 (FIG. 6) involves making an efficient 
determination of the correct table index iythat corresponds to data word fT^for v= 1, 
TV, where the table index /V is a function of start lag At start 93 (FIG. 6), actual PRN 
code chip length 65 (FIG. 6), v, and prompt code chips C(OX C(1X C(2), 
C(A/+1) 91 (FIG. 6). The chip value C(0) 94 (FIG. 6) is needed in order to specify the 
late code 69C (FIG. 2D) at the initial few samples of the accumulation, and the chip 
value C(A/+7) 96 (FIG. 6) is needed to specify the early code 69B (FIG. 2D) at the 
final few samples. Additional constants that can be used in order to determine the /V 
indices are j^s, w^, kmin, kmax, and nominal chip length Atcnom, which has been used 
to generate the three x(i) tables. 

[000S5] The first step of the index calculation procedure pre-computes and 
stores a table of candidate integers for the final summation term that appears on the 
right-hand side of equation (36). This table takes the form: 

Mm)= i;C(//+7-Z-l)2^-'/ for//=l,2,3,...,(M+Z+l) (43) 

7=1 

This computation requires the undefined chip values C(-L+l), C(-L+2), C(-L+3), 
C(-l), and C(M+2), C(M+3X C(M+4), C(M+L). The value zero can be used for 
each of these undefined chips because they can affect the over-sampled codes only for 
the first nexo samples of data word Wi 95 (FIG. 6) or for the last w^u^ samples of data 
word Wn 99 (FIG. 6), none of which are part of the accumulation. The table of 
equation (43) can be constructed by using the following iterative procedure: 
zl/(l) = C(0) (44a) 
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= wa42zl/(//-l),2^] + C(//-l) for// = 2, 3, 4, (A/+2) (44b) 

Ai(ju) = mod[2Ai{ju - 1),2^] for // = (A/+3), (A/+4), (A/+Z+1) (44c) 

Note that the mod(2x ,2^) operation in the latter two equations can be replaced by a 
single truncated leftward bit shift, 

[00086] In many cases prompt code chips 91 (FIG. 6) C(0), C(l), C(2), . . . can 
be generated as the output of a feedback shift register or a system of such registers. 
For example, the new GPS civilian L2 signals can be generated this way. In this case, 
each iteration of equation (44b) can be interleaved with an iteration of the shift register 
calculations. Shift-register generation of PRN codes is well-known in the art. 
[00087] An alternative to building up the previously-described table is to 
calculate the index component only for one data word at a time. Suppose that Aiyis 
the correct index component for data word Wi^ and that /^y is the auxiliary index that 
would have been used to determine Aiy ^om the Ai(ji) table had the table existed. In 
order to calculate Aiyi.\ for data word Wi^u /^vn is computed (procedure defined 
herein), feedback shift register calculations that generate C(//v), C(jUy^lX C(//^H-2), 
Cijiv^i -1) are iterated, and the resulting chip values are used to perform (jUy^i - //v) 
iterations of equations (44b) or (44c). 

[00088] Determination of the correct index into the JCp(/), Xemizm(i)y and Xemt2s(i) 
tables for data word JVy can be reduced to the determination of two quantities. One is 
the time offset index ^i,that causes AtokySrom equation (32) to match the true time 
offset for data word Wy as closely as possible. The other quantity is the auxiliary table 
index juy. It constitutes an index for the sequence of actual code chips that are 
associated with data word Wy. Given these two quantities, the correct index for the 
three x(i) tables is 

= 1 + (ky - k^i^) X 2^ + Aiitiy) for v= 1, 2, 3 N (45) 

[00089] The auxiliary index piy is determined by the position of the Wy data 
word relative to the PRN code chip sequence. Once that position has been 
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ascertained, the index ky can be calculated from the position relative to the samples 
of the L code chips that are associated with the index //v. 

[00090] A time integer can keep track of the number of fine-scale time units in 
a given interval. The fine-scale time unit is a small fraction of the sample interval 63 
(FIG. 2D), At,: 

Atf=^ (46) 
my 

where m/is the integer number of fine-scale time intervals per sample interval 63 
(FIG. 2D), Ats. This number is chosen large enough, for example /w/> 2mN, to 
preclude any significant build-up of timing errors during an accumulation interval due 
to the finite time resolution At/. Nis the number of data words in the accumulation 
interval. The calculation of the values over one accumulation interval involves 
approximately AT iterative time increments, each of which has a resolution of Atf. If m/ 
obeys the inequality given above, then the cumulative timing errors due to the finite 
precision At/ will be less than the timing error caused by the finite timing precision of 
the x(i) tables. Normally it is possible to make m/much larger than 2mN and still 
keep all of the relevant calculations Avithin the size limits of a 32-bit signed integer. If 
m/is a power of two, a rightward bit shift: operation can be used to implement integer 
division by mf Time unit Atf can be used to define an integer that approximately 
keeps track of the code/sample time offset zl/o^^for data word Wy: 

k^^ = round{^) = round{-^^) (47) 

where the roundQ fijnction rounds up or down to the nearest integer. The time lag 67 
(FIG. 2D), zl/ov, is the amount by which the end time of PRN code chip C(jUvL) lags 
the first sample time of data word Wy . The algorithm that iteratively determines k/y 
tries to keep the relationship in equation (47) exact, but using only integer operations 
can allow small errors to build up. Note that kfjmf^ kjm, as implied by a comparison 
of equations (32) and (47). This relationship can be used to determine ky from an 
iteratively determined kfy. Several constants are required by the iterative procedure 



43 



that detornines kfv, kv, and The first five constants are used to account for the 
difference between the nominal chip length Alcnom, used to generate the tables, and 
the actual chip length 65 (FIG. 2D), zl/cused in the accumulation: 

k fmid = round[(n^ - l)m y/2] (48a) 



"■cnom 



Mr 



a 



fixO 



- ceil 



^finid 



m 



sign{X) 



f 1 '^^ ^^c^^^cnom 

f>fix^\ ceil[log{af^/X)/log{2)\ 



^ fix - f'Ound{Xb jij^) 

where the signO function returns +1 if its input argument is positive, zero if the 
argument is zero, and -1 if the argument is negative. The index kfinid is approximately 
half the length of a data word as measured in units of zl//seconds. During an 
accumulation, the rational factor Qfijbfix gets multiplied by the time offset between the 
end time of the first code chip and the midpoint of the data word. The result is a time 
perturbation that removes the average effect of the difference between the actual and 
nominal PRN code chipping rates. The time perturbation can be used to compute a 
corrected kfy value: 



(48b) 
(48c) 

(48d) 
(48e) 



^ Mix ikfy)^kfy^ round 



a 



fix 



^fiX 



(49) 



[00091] Equation (48d) picks bfix to equal a power of two so that the integer 
division by bfix in equation (49) can be accomplished using a rightward bit shift 
operation. The roundQ operation in equation (49) can be accomplished as part of the 
division if one first adds sign{aflx)xbfix/2 to the quantity {kfinid-kf^xafix before 
performing the rightward bit shift that constitutes division by 6/^. This approach can 
give the correct kfypx because the signs of {kjmid-kfy) and bfix are both positive and 
because the rightward bit shift has the effect of rounding the signed division result 
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towards zero. An alternate implementation of the round function could be used for 
applications that do not guarantee kfinid > kfy. Such applications are normally 
associated with L<2 PRN code chips per data word. 

[00092] Five additional constants can be used to define the kf^ and jUy iterations: 



Mjc - round ( 



Akfty^=round( 



kfmin = round 



k f max = round 



^fLtypAtc 
At, 



)'n^mf 




(50a) 
(50b) 
(50c) 

(50d) 
(50e) 



The constant Ltyp is the typical number of code chips per data word. It is the nominal 
increment to //i^per data word. The constant Akfc equals the number of fine-scale time 
intervals per PRN code chip. The constant Akfc is used to adjust kfy up or down if kfy 
falls outside of the limits: kfmm < kfy < kfinax- The constant Akjfyp is the nominal 
increment to kfy per data word. The limits kfinin and kfinax are approximately the limits 
kmin and k^ax fi*om equations (33a) and (33b) re-scaled to the new fine time scale and 
adjusted for the difference between the nominal code chipping rate of the x(i) tables 
and the actual chipping rate of the accumulation. The extra -2 term on the right-hand 
side of equation (33a) is compensated for by the increment to kmm on the right-hand 
side of equation (50d) and the decrement to k^ax on the right-hand side of equation 
(50e), The original -2 term and the increment and decrement have been included 
because they ensure that values which respect the limits in equation (50d) and (50e) 
are transformed into k values that respect the limits in equations (33a) and (33b). 
[00093] The iteration begins by initializing kf\ and fx\ for the first data word. 
The nominal initial values are: 
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kf\nom= round 



At 



start 



Mr 



+ 1 + floor 



-At 



start 



At. 



Ml 



nom 



= floor 



At 



start 



(51a) 
(51b) 



c J 



It is possible that kf\nom from equation (51a) can violate its upper limit kfmax' 
Therefore, the following conditional adjustment can be implemented in order to finish 
the initialization. 



^fXnom 



^ ^fXnom — ^fmax 



Ml =' 



^flnom 
Ml nom 



<ki 



^ ^flnom — ^fmax 



mom 



[Mlnom -1 ^^fmax <^fl 

Given this initialization, the calculation of {^fir,H'i)^ (*y3,//3), {kp^M^:^ ikfN,MN) 
proceeds according to the following iteration: 

fvnom = ^ /(v-1) + ^ ftyp 
Mvnom = /^(v-1) ^typ 



for v=2, 3, 4,..., AT 
for 3, 4,,..,7\^ 



^ fvnom ^fc fvnom ^ ^/min 



fvnom 



\^kr < kr kjr 

'^fmm — ' fvnom — ' /max 



for v=2, 3,4,...,A^ 



Mvnom 
Mvnom 
Mvnom 



+ 1 if A. „ < kf- 

if >tr < kr < kr 

Vmm — fvnom /max 

1 if ^ ^ ^ /".M^ — 

/max jvnom 



for v=2, 3, 4,...,JV 



The table look-up calculations fuiish with the computation of ^v;, 'm, and the actual 
table look-ups: 



ky = round 



m 



f 



for »/= 1, 2, 3 



(52a) 
(52b) 



(53a) 

(53b) 

(54a) 



(54b) 



(55) 



[00094] The roundQ operation in equation (55) can be implemented by adding 
m/2 to mxkftfix(kfy) before the rightward bit shift that constitutes division by w/. The 
result of the division will be the correct value of ky for any sign of k/^k/^ if the 
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computer works with 2's compliment notation for signed integers and if the rightward 
bit shift fills in fi-om the left with the 2*s compliment sign bit, i.e., with the left-most 
bit. 

[00095] Given A:^,from equation (55) and jUyfrom equation (54b), one can use 
equation (45) to compute /V. This value, in turn, can be used to index into the tables to 
determine the Prompt PRN code sign 29A (FIGs. 2D and 3 A), Xpv, the early-minus- 
late PRN code zero mask 35B (FIGs. 2D and 3 A), Xemtmvy and the early-minus-late 
PRN code sign 35 A (FIGs. 2D and 3 A), Xemnsv, that correspond to data word Wy. 
Xpy = Xp{iy) for v= 1, 2, 3 N (56a) 

Xemkmv = Xemhm{Q for V= 1, 2, 3 (56b) 
Xemtlsv = Xemniii^) tOT V= 1, 2, 3 N (56c) 

[00096] The conditionals in equations (54a) and (54b) can be reduced to a 
single conditional per data word during normal operation to improve efficiency. This 
can be done because the sign of Akftyp in equation (53a) is fixed for a given 
accumulation interval (Normally the sign of Akftyp does not vary fi-om accumulation 
interval to accumulation interval or from channel to channel for a given receiver 
because the only variable quantity that affects Akftyp is actual chip length 65 (FIG. 2D), 
Ate, which normally does not vary significantly .) If Akftyp < 0, then the proper formula 
for determining kj^^ and ju^cm be chosen by considering the inequality kfynom < kftnm- 

Conversely, if Akftyp > 0, then the proper formula can be determined by considering the 
inequality kfynom > kfinax^ The decision about which condition to check can be made at 
the beginning of the accumulation because Akftyp is calculated prior to execution of the 
iteration in equations (53a)-(56c). 

[00097] When using a processor that creates instruction pipehnes, "if 
statements can disrupt the pipeline. In this case equations (54a) and (54b) can be 
replaced with the following computations: 



min[0, sign{k y^^^ - k )] Ak^pKO 

0 ^ftyp^^ for v= 2, 3, 4,..., TV (57a) 

max[0, sign(k j^^^ • k )] Akfyp>0 



47 



fv ~'^fvnom 



for v=2,3,4,...,N 



(57b) 



for v=2, 3,4, ...,A^ 



(57c) 



The wmO and moxi^ functions return, respectively, the minimum or maximum of their 
two input arguments. The variable 7/vis normally zero, in which case equations (57b) 
and (57c) leave kfy equal to kfynom and equal to ^vnom^ The value of rifv is -1 ifAk/fyp 
< 0 and kfvnom < kfinin, and +1 if Mftyp > 0 and ^/wjow > kfinax- In both of these cases rjfy 
causes equation (57b) and (57c) to perform the necessary adjustments to kfy and //v. 
Note that eflFicient code may not execute the conditional in equation (57a) once per 
data word. Instead, its accumulation iterations could be performed in one of three 
different iterative loops, depending on the value of Akftyp, Additional economies can 
be had in the first and third conditional clauses of equation (57a). The value of -rffy for 
the first condition is equal to the sign bit of the 2's compliment representation of kfynom 
- kfinin^ Similarly, -^rjfy for the third condition is equal to the sign bit of the 2's 
compliment representation of kfrua - kfynom- In either case, ijfy (or its negative) can be 
computed in two operations. 

[00093] Summarizing real-time over-sampled PRN code generator 3 OA (FIG. 
3 A) and referring now to FIG. 7, to compute prompt PRN code 29 (FIG. 3 A) and 
early-minus-late PRN code 35 (FIG. 3 A) for an entire accumulation interval, the 
method includes the steps of iterating equations (44a)-(44c) (method step 201, FIG. 7) 
to construct the table of Ai(ji) values. The method further includes the step of 
computing the auxiliary constants (method step 203, FIG. 7) in equations (48a)-(48e) 
and (50a)-(50e). The method further includes the step of initializing kf\ and ju\ 
(method step 205, FIG. 7) by evaluating equations (51a)-(52b). The method further 
includes the step of iterating equations (53a), (53b), (57a)-(57c), (49), (55), and (45) 
(method step 207, FIG. 7) to compute, for each iteration, kfynom, Mvnom, ^fv, kfy, ^y, ky, 
and iy. The method fiirther includes the step of iterating equations (56a)-(56c) 
(method step 209, FIG. 7) to compute, for each iteration, Xpy, Xemizmv, and Xemnsv- 
[00099] As mentioned already, it may prove efficient to interleave the equations 
(44a-c) iterations and the accompanying shift register iterations between the iterations 
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that compute kfynom through Xemtisv^ In this scenario //„ can be computed from equation 
(57c). Afterwards, the shift register iterations that generate code chips C(//v^i -1) 
though C{jLiv'\) can be performed, and these chip values can be used to iterate 
equations (44a-c) from /i^^x to fiy xn order to determine /^(jUy) from zl/(//v^i). 
[000100] The software correlator 19 (FIG. 1) of the present invention can 
advantageously be easily modified to work with signals at different frequencies, new 
PRN codes, or even signals for different types of devices. Thus, the same hardware 
could use the software correlator 19 (FIG. 1) to implement such devices as a GPS 
receiver, a cell phone, or both. To allow for new codes, new frequencies, and new 
types of functionality, small changes can be made in the software correlator 19 (FIG. 
1), or different versions of the software correlator 19 (FIG. 1) can be run on the same 
processor. The changes involve using a different baseband mixing frequency and a 
different PRN code in the correlation, and perhaps changes that would provide the 
new signals of interest to the software correlator 19 (FIG. 1). In order for the present 
invention to work with signals at different frequencies, new PRN codes, or signals for 
different devices, two fundamental changes need to be made. First, the baseband 
mixing frequency must be tailored to that of the signal data 21, which also involves 
pre-computing and storing sine and cosine tables at this new frequency. Second, new 
pre-computed over-sampled PRN code tables 28 (FIG. 3 A) must be constructed. The 
size of the new tables should match the over-sampled accumulation period, or at least 
one over-sampled period of the PRN code. As an alternative to generating new pre- 
computed over-sampled PRN code tables 28 (FIG. 3 A), the new PRN codes can be 
generated in real-time by over-sampled PRN code generator 30A (FIG. 3 A). Also, the 
system and method of the present invention could be implemented within systems 
such as GLONASS receivers, cell phones and cell base stations, pagers, wireless 
Ethernet (e.g. 802.11x standards), Bluetooth™, Blackberry® wu-eless internet devices, 
and satellite radio/phones (e.g. INMARSAT®). In fact, the system and method of the 
present invention are applicable to any sort of telecommunication system/device that 
uses spread spectrum, code division multiple access (CDMA) pseudo random number 
codes for the transmission of information, either wired or wireless. 



49 



[000101] Referring now to FIG. 8, navigation calculations require measured 
values of the PRN code phase 55, carrier phase, and carrier frequency. The 
measurements for all tracked satellites must be taken at exactly the same time. A time 
interval counter (TIC) fonction provides a periodic timing scheme to synchronize 
these measurements at time 56. At time tj^^ the TIC Sanction latches all of the 

PRN code phase 55, carrier phases, and carrier frequencies along with the code epoch 
counters, and software correlator 19 (FIG. 1) makes these available to application- 
specific code 15 (FIG. 1), for example, GPS receiver software. GPS receiver software 
uses the code phase and epoch counters to compute the pseudo range to each satellite. 
Software correlator 19 (FIG. 1) keeps track of the code and carrier phase of each 
signal as determined by the code chipping rate and the carrier Doppler shift inputs. 
The quantity the estimated code chipping rate of software receiver 10 for satellite 

j during its k^^ PRN code period, can be determined either by an acquisition search 

procedure, or if tracking, by a delay-locked loop. Likewise, a>Dopmk > the associated 

carrier Doppler shift, can be defined by an acquisition procedure or, if tracking has 
commenced, by a phase-locked loop or a frequency-locked loop. These 
determinations are made by application-specific code 15 (FIG. 1). Software correlator 
1 9 (FIG. 1) can use these two frequencies to update quantities that keep track of its 
code and carrier phases according to the formulas: 



1023 



(58) 



^jk+l " ^jk ^Oomk+l(Tjk+l - Tjk) 



(59) 



[000102] Software correlator 19 (FIG. 1) can keep a running track of these 
quantities and can initialize these iterations as part of the signal acquisition 
calculations that it carries out in conjunction with application-specific code 15 (FIG. 
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1). The quantities tjq and ^josre either sent to software correlator 19 by application- 
specific code 15, or they are initialized arbitrarily by software correlator 19 and 
application-specific code 15 executes feedback control of fcjk ^Doppjk ^ force 

the sequences defined by equations (58) and (59) to converge to appropriate values. 
Information about the previously-described conventional method can be found in 

Dierendonck. 

[000103] The TIC time t^^ 56 (FIG. 8) can occur at, for example, the 
millisecond boundaries of the receiver clock. At each time /^^ 56, the PRN code 
phase 55 (FIG. 8) of each signal is computed in the following manner: 

y^^,=1023 (^^7r) (60) 

Tjk+2 Tjk+\ 

where ij/^^^^ is the PRN code phase 55 (FIG. 8) in chips of signal J at TIC time t^^. 

56 (FIG. 8). The epoch counters, which are simply a running total of the number of 
code periods 57 (FIG. 8), are incremented at each code start/stop time. 
[000104] The carrier phase calculation at time tj^^ 56 (FIG. 8) is similar to the 
PRN code phase 55 (FIG. 8) calculation: 

^jnc " ^ jM ^DomMifnc ~ f y*+i) (6 1 ) 

where ^^^^^ is the carrier phase at time t^^. 56. The Doppler shift that gets returned at 
time t^,c 56 (HG, 8) is a)Doppik^x • 

[000105] With respect to the performance of the system and method of the 
present invention, a sample screen-shot from the illustrative embodiment of the 
present invention is provided in Table 8. This table illustratively shows the tracking of 
nine channels. The roof-mounted LI antenna of the illustrative embodiment can have a 
pre-amp with 26dB of gain. The software correlator 19 (FIG. 1) of the present 
invention can provide positional accuracy on the order of 10-15 meters when working 
in conjunction with application specific software 15 (FIG 1). 
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Table 8. Numerical GPS tracking data output by the illustrative embodiment. 
[000106] Two comparison tests illustrate the performance of the system and 
method of the present invention. In the first test, a first configuration includes a 
MITEL® GP2021 hardware correlator, but is in all other ways identical to a second 
configuration that includes the software correlator 19 (FIG. 1) of the present invention. 
The two configurations differ in SNR by less that IdB and in navigation solutions by 
no more than S-10 meters. In the second test, timing studies using the system of the 
present invention show that processing six channels uses only about 20% of the 
processor's capacity, while Akos 2001a report a real-time software GPS receiver that 
would require 100 % of the capacity a 1.73 GHz microprocessor to implement a 6- 
channel GPS receiver when processing data fi^om an RF fi-ont-end with a sampling 
frequency of 5.714 MHz. 

[000107] Referring now to FIGs. 9A and 9B, among other indicators that could 
assess the accuracy of the PRN code generated by real-time over-sampled PRN code 
generator 30A (FIG. 3A), which includes prompt PRN code sign 29A (FIG. 3 A), 
early-minus-late PRN code sign 35A (FIG. 3A), and early-minus-late PRN code zero 
mask 35B (FIG. 3 A), is the low distortion of the generated codes versus the true 
codes. FIG. 9A, generated for prompt code comparisons, shows juxtaposed plots of 
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I 
I 

I the autocorrelation function 80 (FIG. 9A) of the sampled true code and the cross 

correlation function 82 (FIG. 9A) between the sampled true code and the sampled 
code as generated from the new Xp table. FIG. 9B shows a similar comparison for 
early-minus-late PRN code 35 (FIG. 3A). In either case, the nominal chipping 
frequency is fcnom = l/^tcnom = 1 023 MHz, the sampling frequency is fs = l/Ats = 
12. 199 MHz, and the code timing resolution of the table is /w = 12, i.e., 1/12"^ of a 
sample. This resolution translates into approximately 1/143"^ of a code chip. The 
actual chipping rate differs from the nominal by (^ - fcnom) =15 Hz. This Doppler 
shift of the code chipping rate corresponds to a significant non-zero range rate 
between the receiver and the transmitter, 438 m/sec. The correlations are accumulated 
over 10230 code chips, and the code is a time-multiplexed version of a pair of the new 
GPS civilian L2 CL and CM codes. The data word indexing calculations for this 
example use a code chip start time resolution of AtJlO^, i.e., mf= 10^. 
[000108J It is obvious from FIGs. 9A and 9B that prompt PRN code 29 (FIG. 
3A) and early-minus-late PRN code 35 (FIG. 3A) distortion is very small. The 
correlations produced using the new tables are virtually identical to those produced 
from the exact code. The low distortion of the new approach is best characterized by 
two parameters: the amount by which the correlation peak of the prompt PRN code 
sign 29A (FIG. 3 A) droops below one and the offset of the zero-crossing time of the 
early-minus-late PRN code 35 (FIG. 3 A). The former metric characterizes the power 
loss of the new approach, and the latter metric characterizes the net timing error. The 
droop of the prompt peak is only 0.3% of the nominal amplitude, which translates into 
a 0.03 dB loss. The timing distortion of the new code is less than 4x10"^ code chips. 
This distortion is very small; it translates into about 0. 1 m of GPS range measurement 
error. 

[000109] It may seem paradoxical that the code timing error is only 4x 10"^ code 
chips when the code timing granularity of the x(/) tables is AtJ{mAtc) = 7x10'^ code 
chips. The resolution of this paradox lies in the averaging effect of the accumulations. 
The length of a code chip equals 143,099269 code offset time grid intervals for the 
example shown in FIGs. 9A and 9B. The non-integer nature of this number causes the 
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code offset errors of the x(i) tables to get dithered as the accumulation works its way 
through successive data words. This dithering tends to average out the table 
granularity errors, and this averaging can reduce the net timing error by an order of 
magnitude or more, as shown in FIG. 9B. 

[000110] Although the invention has been described with respect to various 
embodiments, it should be realized this invention is also capable of a wide variety of 
further and other embodiments. 



